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Abstract 

We provide an algorithm for solving the initial value constraints of 
GR using the rectangular parametrization of the CDJ matrix. This 
paper focuses on the Gauss' law constraint and the criteria for its 
integrability. The Gauss' law constraint takes as its input the cho- 
sen configuration, and for an output reduces the CDJ matrix to its 
physical degrees of freedom. This in turn provides an input into the 
Hamiltonian constraint for which we provide a formal solution by ex- 
pansion about the Kodama state in the general case. We establish 
criteria for integrability of the Gauss' law constraint for given starting 
configurations and we provide several examples. 
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1 Introduction 



The canonical formulation of general relativity in the metric representation 
produces a totally constrained system as a consequence of diffeomorphsim 
invariance. The Hamiltonian consists of a linear combination of first class 
constraints = {H, Hi), respectively the Hamiltonian and diffeomorphism 
constraints. These constraints Hf^ have thus far turned out to be intractable 
in the metric representation due to their nonpolynomial structure in the ba- 
sic variables. A major development occured in 1988 with the introduction 
of the Ashtekar variables (see e.g. [1],[2],[3]), which led to a simplification 
of the initial value constraints into polynomial form. The basic phase space 
variables in the Ashtekar formalism are a left-handed SU{2)- Ashtekar con- 
nection and its canonically conjugate variable, a densitized triad 5* . The 
Ashtekar variables effectively enlarge the metric phase space il, essentially 
by embedding it into the phase space of a 50(3) Yang-Mills theory. A rem- 
nant of this embedding is the inclusion of the Gauss' law constraint Ga in 
the list of initial value constraints i?^ — t- (i?^, Ga)- 

A general solution to the initial value constraints problem of general 
relativity entails finding a projection from the unconstrained phase space 
to the constraint surface defining the physical degrees of freedom. It is 
currently an open problem which has been transformed by various authors 
into a form more amenable to physical interpretation [4], [5]. It was shown 
in [6] that the algebraic subset of the initial value constraints H/j^ could be 
solved by a substitution known as the CDJ Ansatz. 

= ^aeBl (1) 

The CDJ matrix takes its values in two copies of SU{2)-, and the 
magnetic field Bl must be nondegenerate. The constraints can in general 
be solved whether for just pure gravity, or to include matter couplings as 
demonstrated by Thiemann [4], by use of the CDJ Ansatz^ 

The main obstacle to the claim of a complete solution to the initial value 
problem, as noted and by the interpretation of the present author, is the 
disparity in the form of the Gauss' law constraint Ga in relation to the dif- 
feomorphism and Hamiltonian constraints i?^ = {H,Hi). Whilst i?^ are 
simple algebraic equations, Ga is a set of differential equations. Neither 
has the existence of solutions to Gauss' law conjunction with the remaining 

^The index convention is that lowercase Latin symbols a,b,c, . . . from the beginning 
of the alphabet will denote internal left-handed 5*17(2) _ indices, and from the middle of 
the alphabet i, j, k will denote spatial indices. We will often denote a position vector in 
three space at fixed time t by r. 
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constraints been addressed in [5]. Thiemann in [4] illustrates an algebraic 
solution of Ga involving an alternate parametrization of the basic variables. 
However, the result is to convert constraint into a differential equation 
in this parametrization, hence presenting the same dilemma albeit in a dif- 
ferent form which would seem to imply a kind of duality between gauge 
transformations and diffcomorphisms. Other approaches to Ga include [7], 
though the Hamiltonian constraint remains unresolved. 

Our proposal for implementing the initial value constraints uses an in- 
stanton representation of Plebanski gravity as a starting point. The initial 
value constraints in the instanton representation are given by 

H = (detS)VVdet*(^ ((tr*)^ - tr*^) + Adet*) = (2) 

for the Hamiltonian constraint, where A is the cosmological constant. The 
diffeomorphism constraint is given by 

= e,jkBiB^,^ae = 0, (3) 

where i?* is the magnetic field derived from the Ashtekar connection A", 
given by 

Bi = e'^^djAl + ^e'^^P^''A]Al. (4) 
The Gauss' law constraint is given by 

Ga=Ve{*ae} + C„^^*/g, (5) 

where Ve = B^di is a triple of vector fields constructed from 5*, and 



= {fabfSge + fbgeSaf)Cbe (6) 

where Cj,e = A^Bl is the so-called magnetic helicity density matrix.^ The 
vector fields Vq satsify the commutator bracket 



[Va,V6] =/i^6Ve, (7) 

where the structure functions arc given by 



^See Paper VI regarding the background and physical interpretation behind the mag- 
netic helicity density matrix. 
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(8) 



The constraints (2), (3) and (5) can be transformed into the Ashtekar 
variables for nondegenerate i?* and nondegenerate using (1). The 
present paper will be concerned mainly with the Gauss' law constraint Ga- 
To obtain the initial value constraints in the Ashtekar variables, one could 
substitute the inverted form of the CDJ Ansatz 



into (2), (3) and (5). This equivalence holds only when the variables are 
nondegenerate. 

The organization of this paper is as follows. In the present section we 
have written the initial value constraints of GR on (^'ae, Af), the phase space 
of the instanton representation. The constraints in and of themselves are 
ultimately constraints on '^ae, the momentum space variables of this repre- 
sentation. Section 2 expresses the Gauss' law constraint in operator-matrix 
form as a triple of differential equations on ^ae- These equations are pre- 
sented in coordinate-invariant form using vector fields v^ = B^di tangent to 
7, the congruence of their integral curves chosen to fill 3-space S. We present 
a procedure which inverts the matrix of these vector fields, taking their non- 
commutativity into account. The result is a propagator from the diagonal 
(anisotropy) to the off-diagonal (shear) elements of ^ae, which constitutes 
a reduction in three degrees of freedom. Section 3 provides arguments for 
the invcrtibility of these vector fields from different perspectives, based 
on the chosen configuration A'^. In section 4 we provide several examples 
for various configurations Af where we invert the matrix of vector fields and 
solve the Gauss' law constraint. Certain configurations are integrable while 
others are not. By integr ability we mean in the sense that the solution is 
independent of the particular path taken along the integral curves 7, and 
depends on the starting and ending points. For the latter nonintegrable con- 
figurations, the particular path must be specified in the solution since the 
solution may be path-dependent. In either case, the algorithm results in a 
reduction in '^ae by three D.O.F. Section 5 treats the inversion of the Gauss' 
law differential operators for the general case without regard for integrabil- 
ity, and section 6 provides a formal solution to the initial value constraints 
by expansion about the solution corresponding to the Kodama state. 

It will be convenient to regard the basic physical degrees of freedom 
which are independent of the gauge and of the coordinate system as resid- 
ing within three elements of ^ae, the diagonal elements. The coordinate- 
dependent, gauge-dependent degrees of freedom are contained within the 



(9) 
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connection Af. But it is also convenient to regard these latter degrees of free- 
dom as more fundamentally residing within 7 which is coordinate-invariant. 
Then the connection and all quantities needed to invert the Gauss' law con- 
straint become derived upon the choice of a coordinate system. 
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2 Gauss' law constraint in rectangular form 

The purpose of the Gauss' law constraint in the instanton representation is, 
starting from a given configuration Af, to estabhsh a map from three CD J 
matrix elements to the remaining six, which constitutes a reduction in 
three degrees of freedom. Starting from the Gauss' law constraint Vei^ae}'^ 
Ci^^ fg = 0, we expand into constituents 



where we have used the Einstein summation convention. Equation (10) is a 
set of three linear first-order differential equations in nine unknowns, which 
implies a reduction of by three degrees of freedom. 

It will be necessary to disentangle the CDJ matrix elements in order 
to assess the nature of the vector spaces being acted upon by the Gauss' 
law differential operators. Expanding (10) to first disentangle the diagonal 
(anisotropy) CDJ matrix elements, we have 



Vl(^ll) + V2(^12) + V3(^i3) + Cl^'^fg = 0; 
Vl(^'2l) + V2(^'22) + V3(^'23) + Cl^"^ fg = 0] 
Vl(*3l) + V2(*32) + V3(*33) + C'/^^/g = 0, 



(10) 





/,fl7^2,2 




/,97^3,3 



The full expansion of (10) reads 




+ *21 + ^^^-31 + Cf *23 + Ct^: 



32 




(12) 



for the first equation, 




+ [C2^'^'l2 + Cf ^32 + C2^'*13 + C|'*3l] + (v3(^'23) + C|'^: 



23 



;) = (13) 
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for the second equation, and 



(vi(^'3l) + *3l) + (v2(*32) + Cf ^'32) + (v3(^'33) + Cf ^33) 



CP^12 + CP^13 + cf *23 + C|l*21 



(14) 



for the third equation. This system of equations can be more conveniently 
put into matrix form 



+ 



+ 



Cf Cf' (V3 + C| ; 

(vi + C7|i) Cf 

cf (V2 + cf ) 

We next expand the matrix equation into diagonal and off-diagonal parts, 
using = ^(oe) + ^[ae]- This yields, for the parts involving vector fields, 








1 


^-22 




^ ^33 / 




f *12 \ 




^23 




V *31 / 




Vl 














V2 




a 










V3 / 









-V3 \ 


/ *[12] 


° 


*[23] 


Vl / 


V *[31] 



2.1 Transformation into Cartesian basis 

It will be useful to transform the Gauss' law constraint into a basis of irre- 
ducible parts of SU{2)- ® S'f7(2)_. Let us now decompose the CDJ matrix 
^ae into orthogonal subspaces consisting of diagonal and off-diagonal parts 
respectively, given explicitly by 



*ae = ^/(e^)ae + ^f{^f)ae + ^f{E^)ae- (15) 

The following identifications have been made 
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^-(12) = ^s; *(23) = ^i; *{3i) = *2; 



*[i2] = V's; *[23] = '^i; *[3i] = V'2; 

*11 = *22 = (^2; *33 = (16) 



Equation (15) in matrix form is given by^ 

(ipi - ^^3 *2 - 1p2 

^3 + V'3 ¥^2 ^1 - V'l 

*2 + ^"2 *1 + tpl (^3 

The basis vectors in this orthogonal decomposition (15) are given by 
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eL = ^afSef, = ^ |e„ed|, (17) 



with no summation convention apphed to the first definition in (17). These 
basis elements have the following matrix representation 




Cael = U U -i ; e„e2 = U U U ; e„e3 







1 
eL = I 1 ; 



ae 















1 

















which satisfy the orthogonality relations 




a,e 
a,e 

{E^\ef)=Y,Ke{efr = ^- (18) 

a,e 

^The diagonal parts <p/ have the interpretation of anisotropy components, while the 
off-diagonal symmetric and antisymmetric components have the interpretation 
respectively of shear and rotation components. 
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with corresponding completeness relations 



E edaeeabf + + E ^aeE'J = SlSf. (19) 

d a 

In the new basis, the Gauss' law constraint can then be written in the 
compact form 

Maf^f + Vaf4^f + Hafipf=0, (20) 

where V, M and H are Gauss' law operators with respect to the shear, 
rotation and anisotropy subspaces, given by 

M„e = ^feW/; F„e = e{gW/; Hae = eie^f. (21) 

The part of the constraint depending purely upon vector fields is given by 



Ve{*ae} 




and the part depending upon the helicity density matrix C^g by^ 



/ Ci3 -(C33 + C22) C12 \ / V'3 \ 

Cl<^^fg=\ -(C1I+C33) W^i] 

V -(C22 + C11) C31 C32 I ) 

(Ci3 — 2C31 C22 — C33 2C21 — C12 \ / ^3 \ 
2C32 — C23 C21 — 2Ci2 C33 — Cii I I ^1 I 
Cii — C22 2Ci3 — C'ai C32 — 2C23 / V ^2 / 

(C32 — C23 —Cz2 C23 \ I '^^ \ 

C3I ^^13 — C3I —Ci3 I I ) 

— C21 C12 C21 — C12 / \ (p3 / 

The vector field contributions are matrices whose elements consist of non- 
commuting differential operators. We will need to invert these matrices in 
order to implement the Gauss' law constraint. 

*See Appendix A for the derivation of the individual terms. 



8 



2.2 Formal inversion in the Cartesian representation 

Prom the diffeomorphism constraint (3) , the antisymmetric part of the CD J 
matrix vanishes tjjd = (dae^ae = 0, and the Gauss' law constraint reduces to 
just the symmetric part 



Hae^e + Mae^e = 0. (22) 

Equation (22) is a hnear relation between the shear elements *e and the 
anisotropy elements ifg. We can reduce by three additional D.O.F. by 
solving for ^'e in terms of ipe, which in essence establishes a map (p ^ ^. 
The constraint can be written in the general form 



(Vae + a„e)*e = Qa[0] ^ ^e, (23) 

where wc have defined the operator matrix of vector fields on the shear 
subspace 




with the helicity density matrix insertions 



Cl3 - 


-2C31 


C22 - 


- C'ss 


2C21 


~ C12 


2C32 


— ^23 


C2 I - 


-2Ci2 


C33 - 


-Cu 


Cii 


- C22 


2Cu 


-C31 


C32 - 


2C23 













J 



and on the anisotropy subspace by 



Vi — C[23] — C32 C23 
Qa[0] = - \ C31 V2 - C[3i] -Ci3 

— C21 C12 V3 — C[i2] 

where C^ab] = Cab — C^a- The formal solution to (23) is given by 



(Vea + aea)~^Qa[<^, (24) 

which may be written as an infinite operator expansion, 



Jeeo 



00 



a, 



6364 



n=0 



(25) 
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where v~g acts to the right. We have defined the Gauss' law propagator by 
the operator Jj = — {M^^H)'j, where must be suitably defined. First 

let us carry out the inversion of the part involving just vector fields. 

We would like to express the shear elements ^ as the image of the 
anisotropy elements (p under the map J^. If one chooses v^e to consist 
entirely of vector fields with the helicity density contained in aae, then in 
order to solve Ga one must invert the matrix Vqe- More precisely, we must 
find an operator-valued matrix Uae such that UaeVg/ = (^e/ is the identity 
operator. Since Vae consists of noncommuting vector fields Va, then care 
must be taken with operator ordering when carrying out the inversion. Let 
the required matrix be given by 



Uae = Saiae + 5a2he + ^aSCe, (26) 

where Og, 6e and Cg are operators which remain to be determined. The 
resulting operator product UagVe/ is given by 

oi 02 as \ / V2 V3 \ / aiV2 + 02 Vl 02 V3 + 03 V2 aiV3 + 03 Vi 
hi 62 ^3 Vi V3 = felV2 + h2^l 62 V3 + 63 V2 61 V3 + 63V1 

Cl C2 C3 / \ V2 Vi / \ C1V2 + C2V1 C2V3 + C3V2 C1V3 + C3V1 

We must now choose (og, h^-, Cg) to diagonalize the matrix product, preserving 
the operator ordering. 

First we require that the off-diagonal elements vanish. Starting with the 
first row of UaeVg/, this leads to the equations 

oivs -I- asvi = — > as = -oiVsv^^; 
02 V3 -I- a3V2 = — > a2 = -03V2V3 ^ = oiVsvf ^V2V3 ^. (27) 

The requirement that the off-diagonal elements of Uae'Vef in the second row 
vanish is given by 

61 V2 -I- 62V1 = — > hi = -62V1V2 ^; 
&1V3 -I- 63V1 = — > 63 = -6iV3vJ"^ = 62V1V2 Vsv^^ (28) 

The requirement that the off-diagonal elements of Uae^ef in the third row 
vanish is given by 

C2V3 -I- C3V2 = — > C2 = -C3V2V3 ^; 

C1V2 -I- C2V1 = > Cl = -C2V1V2 ^ = C3V2Vj^ViV2 ^. (29) 
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Having required all off-diagonal parts of Uae^ef to vanish we can can com- 
pute the diagonal parts using (27), (28) and (29). These are given by 



aiV2 -I- a2Vi = ai (v2 -I- vavf V2V3 Vi) ; 
62 V3 -I- 63 V2 = 62 (va -I- viv^V3v5"V2); 
C1V3 -I- C3V1 = C3 (vi -I- V2V3 V1V2 V3) . (30) 

Hence, three elements of Uae, namely oi, 62 and Cg remain undetermined. 
However, the matrix product of UaeVe/ can now be written in the form 



1 V3Vj;V2V3^ -VsV^"^ 

-ViV^^ 1 V1V2 VsV^"^ 

V2V3^ViV2^ — V2V3^ 1 



V2 + V3Vj^ V2V3 ^Vi 

V3 -I- V1V2 ^V3Vj"^V2 

Vi -I- V2V3 ^ViV^^V3 

The undetermined operators appear to the left on both sides. Assuming 
invertibility of these operators we act with Diag{ai^ ,b2^ ,c^^) on both sides 
to obtain the following decomposition 

(V2 V3 ^ ~^ 

VI V3 1 = (l + 7r)-^t;-^M, 

V2 Vi 

where we have defined 

1 V3VfV2V3^ -V3Vf^ 

M = I — ViV^^ 1 V1V2 ^V3V]^^ 




and 



V2V3 V1V2 ^ -V2V3 ^ 



TT = 





as well as the cyclic operators 



7r2 = V2 V3V1 V2V3 Vi; TTs = V3 V1V2 V3V1 V2; TTl = V2V3 V1V2 V3. 



(31) 
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One could then in principle evaluate (25) explicitly term by term, starting 
from a given configuration Af. 

Solving the Gauss' law constraint should reduce by three degrees of 
freedom to the anisotropy subspace iff, given by 



*ae = (eL - ELjf[A])^g = T^^fg. (32) 

The physical interpretation is that upon projection to the Gauss' law con- 
straint surface, the anisotropy basis elements e{e become distorted into T/g 
based upon the configuration Af, in matrix form 



tL= -n -Jl ; Tl= -Jl 1 -J2 .Tl={ -Ji 



1 


-Ji 


-Jl 







-Jl 


Jl 


-Jl 











-Jl 




1 


-Jl 




-Jl 






-Ji 

The shear elements are eliminated and (pf are now expressed in a 'curved' 
basis Tie = Tie[A] which encodes the kinematic effects of the Gauss' law 
constraint for each configuration Af.^ 

2.3 Off-diagonal helicity configurations 

Another possibility to carry out the inversion is to define the insertion matrix 

"ae by 

/ C22-C33 

aae = C33 - Cu 

\ Cu - C22 
which leaves the remaining matrix containing vector fields 



V2 + Ci3 - 2C31 va + 2C21 - C12 

Wae = I Vi + 2C32 - C23 V3 + C21 - 2Ci2 

V2 + 2Ci3 — Csi vi + C32 — 2C23 

The solution to the Gauss' law constraint then amounts to carrying out the 
analogue of (25), namely 



-'eeo 



oo 
n=0 



-1 



-1 



'eoei"eie2 *'^e2e3"e3e4 



lien) ^ena\^afJJf 9^9 = Je^9 



"'\t will be convenient to regard the set of all nonortliogonal charts of E as being the 
fundamental structure. Then the magnetic field B], becomes derived from the vector fields 
tangent to the three linearly independent coordinate directions defining integral curves 
through each point. Hence, there is a correspondence from T/^ to each chart, where the 
transformation from one chart to another induces a passive diffeomorphism. 
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which in the general case may be more manageable. For off-diagonal helicity 
configurations, defined as configurations for Af for which Cn = C22 = C33, 
the inversion can be written in closed form since = 0, which reduces 
equation (33) to just the first term w~g^. One then explicitly writes the 
Gauss' law propagator as J| = wJ^Hae, where 

W2 W3 
= I wi W3 1 = (1 + 7r)-^MJ-^M, 



W2 Wi 



where we have defined 



W3W2 ^W2W3 ^ — WsW^ ^ 



M = I — W1W2 ^ 1 W1W2 ^waW]^ ^ 

W2W3^WiW^''^ — W2W3^ 1 



and 



W2 
W = I W3 

wi 
as well as the cyclic operators^ 



7^2 






















TTl 



7i"2 = wawj" W2W3 wi; 

TTa = W3 ''^WiWj''^W3W]"^W2; 

TTl = W5"^W2W3 ''"WiW^''"W3. (34) 

One has then reduced the problem to the inversion of individual first order 
differential operators of the form w = v + c, where c is a c-number. In 
this case one may be able to find an integrating factor to construct w~^. 
The series (25) formally solves the Gauss' law constraint by construction. 
However, the solution does not exist until it has shown to be well-defined. 
We will put in place the ingredients necessary to properly define all terms 
and their constituent operators, as well as the action within their domains 
of definition, and provide some examples. 



The physical interpretation is best motivated by analogy to the case Va = da, where 
the vector fields are exact partial differentials. Then the vector fields form the structure 
of an abelian group where differentiation are group elements with inverses ~ f dx"", 
which corresponds to autidiffcrcntiation with vanishing boundary data. The cyclic opera- 
tors correspond to making one complete closed path along the edges of a hexagon in group 
space. For the aformentioned case the cyclic operators would act as the identity element 
on the space of functions. In the general case the cyclic operators are not the identity, 
owing to the noncommutativity of the vector fields. 
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3 Hamiltonian flow and integral curves 



A necessary condition for the well-definedness of the expansion (25) hinges 
upon the invertibihty of the vector fields Vq = -B*5i = Xa + N^, where we 
have defined 



Xa = e'^\djAt)di- Na = \e'^^faUA)Al)di = {deiA){A-Xdi (35) 

respectively as the abelian and nonabelian contributions for det^ ^ 0. The 
abelian part of is given by 



Xa = 2 



/dAf _ dA^\ d /dA^ _ dAU d /dA^ _ a^f \ d i 
\ dv dx J dz \ dz dv J dx \ dx dz J dvi ' 



dy dx J dz \ dz dy J dx \ dx dz ) dy 
Rearanging the terms of (36), one obtains 



Xa = Xai + Xa^ + Xa^ (37) 

where we have defined 



^ _ /dA1\ d /dA1\ d 
1 \ dv ^ dz \ dz J dv 



dy ' dz \ dz J dy^ 
dA^\ d /dA^ d 



^'^ \ dz J dx V dx J dz ' 

XAa = (^]--(^)—. (38) 

^ \ dx / dy \ dy / dx 

Each term of (37) exhibits the same form as a Hamiltonian vector field Xjj 
from classical mechanics, given by 



Xh = —--—-. (39) 
dp dq dq dp' 

If one defines a symplectic manifold on phase space Q = {p, q) , with sym- 
plectic two form oj = dpAdq, then the corresponding Hamiltonian equations 
of motion would be given by 



. SH 



'^We have chosen x, y and z as labels for the coordinates, which needn't be Cartesian. 
These symbols can be used for arbitrary 3-dimensional coordinates. 
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The Hamiltonian function H generates evolution in parameter t along a 
curve in $7, and Q is foliated by the orbits of the Hamiltonian flow such that 
a unique curve passes through each point. 

Hence, identify with the components of the connection A" nine Hamil- 
tonian functions A", and ylg on respective phase spaces (y, 2;), {z,x) and 
{x,y). These imply Hamilton's equations of motion 



^ = {^'.^?}; ^ = (x'.^?}, (41) 

with cyclic permutations over i,j,k. Hence, for each Hamiltonian function 
A^ there exists a parameter generating flow along a curve in the two 
dimensional (x-' , x^) hyperplanes dual to each corresponding i direction. We 
would like to envision the Hamiltonian flows as the projection of integral 
curves of Xa onto these hyperplanes. To see that these curves are well- 
defined, re-arrange Xa into the following form 



i=l 

The three dimensional curves (f) for each i, j, k imply the relations 



dAf _ dx^_ a4f _ _dx^ 

'dxT~~dif' 'dx^~~'dtf' ^ ' 



The total variation of Af in 3-space is given by 



dAf = ^dx' + ^dx^ + ^dx^ = ^dx\ (44) 
ax* dx^ dx'^ ax' 

Hence when restricted to the curve parametrized by tf, then Af is constant, 
which is consistent with the 2-dimensional projections implied by (41). The 
vector field Xa can be seen as the composition of all three curves 



d d d _ d 

If the individual contributions could be exponentiated, then (45) has the 
following action on an arbitrary function / 



(Va) 



Abelian 



{/} 



d tx 
dt^^' 



A a tx A a tx , 

^1 e 2 e 



IX 



(46) 
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which ascribes to all three vector fields a flow labelled by the same parameter 
t. The order of the exponentials in (46) does not matter, since they commute 
for t = 0. Hence (45) illustrates the fact that the vector fields form a 
vector space. Since this gives a well-defined fiow, then the inverse of the 
abelian contribution {'v~^) Abelian = exists and admits an interpretation 
of integration along the integral curves (f). 

3.1 Differential geometry of differential forms 

The integrability of the abelian part of Vq can be expressed more succinctly 
using differential forms. This corresponds to three vector potential one- 
forms A" = Afdx^, with field strength given by 

= dA"" + ]^r^''A^ A A". (47) 

The vector field Xa from (35) can be derived from the two form dA"" in the 
following manner. Define r/, the top three form on 3-space S by 

•n = ^Qmndx^ A dx™ A dx". (48) 
The interior product of (48) with Xa is given by 

ixAv) = {V, Xa) = U^hi^ndiA';{dk,dx^ A dx"^ A dx") 
= \e'^^eimndiA''^ (^dx"* A dx^ + 5^dx" A dx^ + S^dx^ A dx'^) . (49) 

All three terms on the right hand side of (49) are equal, summing to the 
abelian part of the field strength, which leads to the relation 

ix^ (r/) = dA'^ d{ix^ (r/)) = d^A^ = 0. (50) 

Therefore, if Xa is integrable, then this must be related to the fact that 
d^A"- = 0. But the vector fields \a = Xa + Na also have a contribution 
from the nonabelian part Na- Therefore we can rephrase the criterion for 
integrability of Vq as a condition on Na, namely that 

d{iNAv))=d{ixAv))=0. (51) 
Using (47) this leads to a sufficient condition for integrability of Vq that 
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/"'"^A^ A dA" = -dif^^A^ A = -^r^^A^ hA'' = d^" (52) 

locally for some one form = ^fdx^. 

There is a wide class of configurations A"' for which this is the case. For 
example, for potentials of the form A'^ = d(p"' for 5(9(3, C)-valued scalars (pa 
we have that Xa = 0, and Na in component form is given by 

v„ = iV„ = ir^^e^^'=(5,r)(9fc/)5,. (53) 

By the above results on differential forms \a is integrable since d(/"'"^^'' A 
d^'^) = 0, which means that exists. Moreover, the helicity density 
matrix 

Cae = AfBi = ^r''^e'^''{djcl>'){dkcl>'^){dir) = 6ae{detd$) (54) 

is diagonal. Hence the c-number insertion matrix aae vanishes and only the 
first term of (25), v~g^, survives. The result is that for the configuration 
Af = di<^°', where = {(f)^ , (fy^ , (f)^) are arbitrary functions in S, the Gauss' 
law constraint should have an explicit solution. 

Equation (51) has the following interpretation, from Cartan's formula 
for differential forms [9] 

Lxari = iXa{dri) + d{ixaV)- (55) 

Using the identity drj = 0, this states that the Lie derivative of the top 
3-form on a 3-dimensional space along the direction of Xa must vanish. 

In the more general case where l\ Ap- ^ d{^f°^<p'^) is not exact, one 
should still attempt to verify the integrability of Ga- The issue of integra- 
bility in relation to the most general magnetic field can best be understood 
from the vantage point of the Bianchi identity, 

dF" + r^^J^ AF^ = 0, (56) 

which wc posit is directly related to the vector fields Vq. The following 
relation ensues, in direct analogy to (50) 

<iri) = F^ (57) 
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where is the total nonabehan field strength. By the Cartan equation we 
have 

L^Ari) = di^M+ivAdv)- (58) 
The second term of (58) vanishes and we are left with 

Lv„ (ri) = dF" = -f^-^A^ A dA" (59) 

where we have used (56) in conjunction with d^A"- = 0. This can be rewritten 
as 

L^AV) = -ir''Cb,)7i- (60) 

Equation (60) implies that a nonvanishing antisymetric part of the helicity 
density could possibly present an obstruction to the integrability of Ga- 

3.2 The various possibilities 

Hence the possibilities can be delineated as follows. If A°- is closed, then 
dA°- = 0, which leads to two subcases according to whether A"- is exact or 
not exact. 

(i) Case(i): A"' closed and exact. In this case one can write A"- = dcj)"- 
for three arbitrary functions 0", and one has that = 0. Hence the only 
remaining issue regards Na, which can be written 

= p^d(j)^ A #^ = d{f"-^''(j)^d(jf) (61) 
which is exact. For this case one has 

v„ = r%fc9,(/5fc(/>^)ai, (62) 

which is integrable. This case for each A"" defines a map between from 
coordinate system {x,y,z) to {(p^ , cp"^ , (p^) , which defnes an equivalence class 
of passive diffeomorphisms. 

Case(ii): Af closed and not exact. In this case Xa is still integrable even 
if not explictly so in terms of global coordinates. It is not possible to judge 
the integrability of Na based on this, since f"-^'^A'^ A dA"^ 7^ in general. So 
one must examine a less restrictive notion of integrability, 

(^fabc^b ^ ^ dir^^Af A ) = (63) 
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with no summation over a. Equation (63) is valid since it corresponds to a 
5-form in three dimensional space S, which must vanish. The result is that 
for A"' closed, the Gauss' law constraint should still be integrable. 

Case(iii): Af not closed {dA'^ / 0). Still, dA'^ is integrable since as 
we have demonstrated it induces a Hamiltonian flow on 2-dimensional dual 
hyperplanes. Intcgrability then hinges upon whether Na is integrable, which 
leads to two subcases. First, if f"-^'^A'' A A"^ is closed then this implies that 



^(^fabc^b ^ ^c) ^ jabc^b ^ ^^c ^ _(^fabc^^-^ ^ g, (64) 

or the antisymmetric part of the magnetic helicity C^ae] ™ist vanish. Cer- 
tainly, as we will demonstrate, there are a wide class of configurations Af 
for which this is the case. 

The second possibility is when Cj^c] 7^ 0, namely that the antisymmetric 
part of the helicity is nonvanishing with Af not closed. While Na does not 
correspond to a closed two form, one can demonstrate integr ability using 
the same criterion as in (63) 



^^abcj^b ^ ^c) ^ ^(^fafg^f ^ ^ q, (65) 

with no summation over a, which is due to the vanishing of a 5-form in three 
dimensional space S. The result is that for generic configurations Af, the 
vector fields comprising the Gauss' law constraint should be integrable. 

Case (iv): We have considered where Xa = and Na ^ 0, which 

yields an infinite set of integrable configurations labelled by three arbit rary 
functions cp"'. Let us now consider the other extreme, where Na = and 
Xa 7^ 0. The nonabelian contributions to are given in component form 

by 



Ni = 


{^4 


- AiAl)d, + {AiAl 


N2 = 




-AlAl)dy + {AlAl 


N3 = 


: {A\Al 


- A\Al)d, + {A\Al 



The total contribution is given by 



-AlAl)dy + {AlAl-AlAl)d,- 

— A^A\)dz + (^2^3 — ^43^2)^^; 

- 4^2)9. + {A\A\ A\Al)dy. (66) 



v„ = X„ + ^cof(^?)^. (67) 

i 

One possible obstruction to intcgrability of the Gauss' law constraint is due 
to Na = coi{Af)di, since we have not shown this contribution to correspond 
to a flow orbit. A sufficient condition for the invertibility of Va is that 



19 



Na = Va. One sees that it is sufficient to require that two rows of cofactors 
to vanish, which ensures that the third row also vanishes. In other words, 
for example taking 

AlAl - AlAl = AlAl - AlAl = AjA^ - AlAf = 0, (68) 

and 

^3^1 ~ ^1^3 — ^1^2 ~ ^2^1 ~ ^2^3 ~ ^3^2 ~ ^ (69) 

yields that = 0, with analogous results for the different combinations 
of rows. The result is that the unconstrained components of Af form a 
five dimensional manifold, with 9!/5!4! = 126 possibilities. We provide one 
particular choice for illustrative purposes 

A^ = ( A\Al{A^-'^ ^ AlAl{Al)-A. 

\ ^3"^l(^3) ^ ^3^2(^3) ^ ^3 / 

If wc require that Na = then Vq is invertible, since we have already shown 
that Xa arc integrable. The components of the magnetic field are given by 

Bl = d2Al - d^Al + AlAl AlAl; 
Bl = d^Al - diAl + AlA\ - A\A\; 
Bl = diAl - d2A\ + A\A\ - A\A\ (70) 

for the diagonal components, 

B\ = d2Al - 83 Al + AlAl - AlAl; 
Bl = dsAl - diAl + AlAj - AiAl; 
Bl = diAl - d2A\ + AlAl - AlAl (71) 

and 

Bf = dsAl - diAl + AlAl - AlAl; 
Bl = diAl - diAl + AlAl - AlA\; 
Bl = d2Al - dsAl + AlAl - AlAl (72) 

for the off-diagonal components. Such configurations, where (detA) = 
with each cofactor identically vanishing, should correspond to a wide class 
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of integrable solutions. Note, with the cofactors vanishing, that the magentic 
heHcity density matrix Cae reduces to 

Cae = e'^'AfdjAl (73) 

Another way to ensure integrabihty of Vq with nonvanishing cofactors is to 
require Na to be proportional to Xa to within a numerical constant. Hence, 
since Xa is already integrable in accordance with the Hamiltonian flow, then 
Na would as well be integrable. 

3.3 Integrability in relation to Bianchi groups 

In the Bianchi models one chooses a connection of the form A'^ = a,aeX^i 
where Oae = £iae(0 is spatially constant, and all spatial dependence is ab- 
sorbed into invariant one forms which satisfy the Maurer-Cartan equation 

dx' = -\c%xf f\x'- (74) 

In (74), are the structure constants for a group of isometries of S, which 
can be decomposed into the form 

C% = rf'^edfg + 5yag^, (75) 

where Og = is the trace of the structure constants. 

The curvature two form for the connection A'^ is given by 

F- = dA'' + ^/«^^^^ A = aaedx" + ^/"''ob/acflX-^ A x' ■ (76) 
For nondegenerate Qae, (76) can be written as 

F« = aaedx' + {deia){a-^r''edfgX^ A x' ■ (77) 
Using (75), (74) can be written 

dx' = —n'^hX^ A - A x'ag. (78) 
This leads to the relation 

edfgX^ AX' = 2{n-y'{-dx' + a^x' A x")- (79) 
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Using the fact that the S'C/(2)_ structure constants for the Ashtekar vari- 
ables are numerically the same as the three dimensional epsilon symbol 
fabc = ^abc, we have 



F" = Uaedx^ + (deta)(a-^)'"' -2{n-'f^dx^ + 2{n-'f''agX^ A x 

= d((a„e - 2(deta)(a-i)«'^(n-i)'^^)x^) + 2(deta)(a-i)«<^(n-i)'^^a,x^ A x^(80) 

The first term on the right hand side of (80) is exact, therefore it is integrable 
and corresponds to a well-defined Hamiltonian flow. The second term is 
not manifestly exact unles = d(j)"' for some which means that the 
Gauss' law constraint might not be solvable if it is not exact. Note for 
ttg = the offending term vanishes, categorizing the Bianchi A models. 
This observation might perhaps shed some light on why Bianchi class A 
minisuperspace models are integrable while class B models, where Qg ^ 
are not. 

An example depicting integrability occurs when Xa and Na are directly 
proportional to one another, namely Na = kXa for some numerical constant 
k. We have effectively analyzed the A; = and A; = oo cases, but for < 
\k\ < oo we have 



dA" = -k^'^'^A^ A A". (81) 

Then Cae = <^ae(A; + l)(det^) which is diagonal. Hence Gauss' law constraint 
is exactly integrable in this case since v~g^ exists and Uae = 0.^ 



3.4 Clebsch decomposition of the potential 

As noted in [13] every three dimensional smooth vector field v may be locally 
decomposed in a region free of singularities as i; = fVg + V/i, and the local 
helicity of v is given by 



H = v- {V xv) = Vf ■ {Vg xVh). (82) 

According to the Frobenius theorem, if the helicity H oi v vanishes every- 
where, then Clebsch functions / and g may be found such that v = fVg. 
Such zero helicity fields are known to be completely integrable in the sense 
of Frobenius. We will now apply the Clebsch decomposition to the non- 
abelian gauge field Af, thought of as a collection of three 3- vectors arranged 

^Equation (81) is essentially the Cartan-Maurer equation where Af are viewed as the 
components of an invariant basis of two forms. This forms essentially the basis for Bianchi 
models in quantum cosmology. 
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in the columns of a three by three matrix. Each column has a Clebsch 
decomposition given by 

= fidgi + dhi; = f2dg2 + dh2\ A^ = f^dgs + dhs, (83) 

namely A" = fadga + dha with no summation over a, where /„, ga and ha 
constitue nine arbitrary smooth functions. For /„ = ga = 0, and ha 7^ 0, 
we have shown already^ that the Gauss' law constraint is integrable. Let us 
now analyze the opposite extreme, namely where ha = with and ga not 
necessarily vanishing. Hence this reduces the Pfaff dimension of Af from 
nine to six.^° For A"- = fadga, the magnetic field is given by 

S» = dfa A dga + ^r'^Mcdgb A dgc. (84) 

The magnetic helicity density matrix is given, written out in terms of com- 
ponents, by 

Cae = e'^''fa{^^ga){^Jfe){^kge) + ^6^^' V^'^a A/c(5^5a) (S,5fe) (5fc<?c) 

= le'''fa{djfe)[{diga){dkge) ' {dige){dkga)] 

+Saeflf2f3{det{dg) = SaeU + (85) 

for some and n constituting respectively, the abelian and the nonabelian 
contributions. One sees that for ha = 0, the diagonal contributions to Cae 
are all equal. The implication is that the matrix ct^e = 0, and hence (33) 
reduces to the first term w~g^, which can in principle be carried out exactly 
in the event that an integrating factor is found. 

This concept can be related to Papers V and VI, starting from an arbi- 
trary potential one form A° = Afdx^. While the curvature of A" is given 

by 



= dA'' + -r^^A" A 0, (86) 

which is nonvanishing for nondegenerate Ashtekar magnetic fields B^, one 
can always find a geometric system for which A^ arc the full-theory analogue 
of Bianchi invariant one forms. Hence there always exist structure functions 

^Recall the results of A^ = dicf)"', whence one makes the identification ha = 4'°' ■ This 
case corresponds to a reduction of the Pfaff dimension of j4" from nine to three. 

^°This will provide as a physical interpretation the analogue of helicity for nonabelian 
gauge fields. 
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dA" + ^K^A'' AA^ = 0. (87) 



Note that one may simply define 



htc = iu-'riM<]h (88) 

where = {A~^y^ for nondegenerate Af. Writing the structure functions 
as a decomposition 



hl = ebcfnf' + ^{Stac-S',a,) (89) 

in direct analogy to minisuperspace, one finds that the components of the 
magnetic helicity density are given by 



Caed^x = A^AdA^ = - (det A) {jf^ + e^^'^a^ . (90) 
3.5 Sufficient condition for integrability 

We will now examine some of the restrictions which can conceivably be 
put on the magnetic field as a consequence of integrability. Specifically, we 
will examine some conditions under which there exist coordinates where the 
vector fields have manifest integrability. The first condition is that 



(h^_dy__dz_ 

for a = 1,2,3. We can write (91) as an exterior differential system 



el = Bldx-Bldy; el = Bldy-Bldz; 9^ = B^dz - B^dx. (92) 

A necessary condition that the one forms in (92) correspond to a reversible 
process is that the topological torsion vanish. For each a we must require 
that 



el A del = [Bl^ - Bl^) dxAdyA dz; 

/ BB^ dB^\ 
el A del = [Bl^ - Bl-gf)dx Ady A dz; 

Ol A del = {bI^ - Bl^)dx A dy A dz. 
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The requirement that each term of (93) vanish leads to the condition 

=a;'"(B|)^ =&'K4) = »- 

Equation (94) leads to the conditions 

B^ B'^ B^ 

-^ = Ca{x,y); -^ = aaix,y); -^ = ba{z,x) (95) 

-"a -"a -"o 

for each a for some functions (oa, ba,Ca)- Additionally, we have that 

aabaCa = 1 Vo. (96) 

We can use (95) to solve for all off-diagonal elements in terms of the diagonal 
elements, obtaining 

/ 1 C2{x,y) b^\x,y) \ / Blix,y,z) 

Bl=\ c^\x,y) 1 a3{y,z) Bi{x,y,z) 

V bi{z,x) a^\y,z) 1 J \ 5|(x,y, 

We can use (96) to eliminate ci and C2, making the re-definition a^^ — > 02 
and b^^ — > 63, obtaining 

/ 1 a2{x,y)b2{z,x) b^{x,y) \ / Sf(x,y,z) 

Bi=\ ai{y,z)b^\z,x) 1 a3(?/,z) Bl{x,y, 

\ bi{z,x) a2{y,z) 1 / V 

We have obtained a decomposition of the magnetic field into a part which 
is restricted by integrability requirements, and a part which is completely 
unrestricted with three degrees of freedom per point. One interesting ques- 
tion is whether this restriction makes sense when plays the dual role of a 
configuration space dynamical variable as in GR. According to the results of 
Paper VI there appears to be no obstruction to any arbitrary GL(3)-valued 
matrix 5^ playing the role of a magnetic field. This is because one can 
always find a connection A*^ such that -B* is the magnetic field, and then the 
helicity density Cae is a derived quantity. So certainly, one can in principle 
find that yields B^ of the above form. 
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4 A few examples 



The Hamiltonian and diffeomorphism constraints can be solved algebraically 
without reference to any coordinate system, and reduce by four degrees 
of freedom upon solution. The Gauss' law contains spatial gradients di 
intermingled with _B* in the form of vector fields V(j. The solution of this 
constraint requires inversion of the differential operators, which requires 
the specification of a configuration of Af under which the inversion will be 
performed. More fundamentally, one must specify which congruence 7 of 
integral curves filling 3-space E one wishes to consider. Note that 7 and 
the tangent vector fields Va are coordinate- independent. The specification 
of a configuration Af then amounts to a choice of coordinates as follows. 
Given 7 and X^, then v^jx*} = determines the components of Vq in a 
particular coordinate system. The requirement that i?* be a magnetic field 
then fixes A^ as one member of an equivalence class of potentials defined 
by such that the Bianchi identity is satisfied. Hence Af in a sense is 
intertwined with the choice of coordinate system. To solve the Gauss' law 
constraint then amounts to the imposition of integrability on , derived 
from 7, which are regarded as the fundamental structure. The CDJ matrix 
*ae suffers a reduction in three D.O.F. upon implementation of Ga, leaving 
only three D.O.F. which can be considered physical in the sense that they 
are independent of the choice of coordinates. 
The Gauss' law constraint can be written as 

^in^fg} = ('^"^(V, - fgbeCbe) + fabfCbg)^fg, (97) 

where the vector fields Vg satisfy (7). The homogeneous reduction can best 
be understood by writing the Ashtekar magnetic field in the form 

Bi = e^^%Al + {detA){A-J,, (98) 

depicting an 'abelian' part involving spatial gradients, plus a nonabelian 
correction. This relation follows from the fact that the structure constants 
fabc = ^abc ^oi SU{2)- are numerically the same as the epsilon symbol ej^fe. 
Likewise, the following suggestive form for the magnetic helicity density 
matrix 

Cbe = e'^'^A^djAl + (detA) (99) 
can be written. Inserting (98) and (99) into (97), we obtain 

w/f = idetA){S-f{A-%d^ + fagf) + via (100) 
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where the remainder , which contains the gradient terms, is given by 

rfj = e'^' [S'^HdjADd, + {Safftge + Segfabf)AtdjAt) . (101) 

4.1 Spatially homogeneous sector of connection superspace 

We now reduce the constraints to the spatially homogeneous sector of con- 
nection superspace, by choosing a configuration for which are spatially 
constant. Note from (99) that the only nonvanishing parts of C^e are 
Cbe = SbeidetA). Setting the spatial gradients to zero, we have rl^ = 
which leads to 

= {det A) {S'^f{A-y^d, + fagf). (102) 
Since fg is symmetric, then the Gauss' law constraint is given by 

w^n^/p} = {dctA){A-')ld,{^ag} = 0. (103) 

For spatially constant S*, we have the relation i3* = (det^)(yl~"^)^. We may 
now rescale the constraint by a factor of det^ ^ 0, since it is now invari- 
ant under rescaling. Making the identification = {A~^yg^di for spatially 
constant Af, we have that [ea,e6] = 0, where 

A; e-i = J dt\ (104) 

Hence the structure functions ^^j, in (8) all vanish. The following relations 
additionally ensue^^ 

[e„, Bb] = [e-\el^] = [e-\eb] = 0. (105) 

Since the matrix elements of the Gauss' law kinetic operator commute, then 
its inversion is analogous to the inversion of a matrix of c-numbers. 
The integral curves of are given by 

^^Minisuperspace, by our definition, means that spatial gradients of all quantities are 
vanishing. By this definition, the spatially homogeneous sector of Ai does not necessarily 

imply minisupcrspacc, since as wc will sec that "tae needn't be spatially constant. In the 
spatially homogeneous sector the identity det_B = {ActAY follows from (98) 

^^We set the constants of antidifTerentiation to zero, so that ^ is the unique inverse 
of Ba- We will adopt the same convention in general when solving the constraint. 
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— = {detA){A-X ^x^ = idetA)iA-yX, (106) 

which are straight hnes in 3-space S. Let us define a 'densitized' version 
_ ^a^det^), of the coordinates. Then the following transformation defines 
a nonorthogonal coordinate system 



= Alx + Aly + Alz; 
? = Afx + Aly + A^z; 

= Alx + Aly + A^z, (107) 

which allows the CD J matrix fg = 'I'jg(P,?,?) to be expressed directly 
in this coordinate system without reference to {x,y, z). 

The matrix form of the Gauss' law constraint (103) is given by 






The solution is given by 



0. 




ei 


















II 










63 / 




{f3 J 



In the spatially homogeneous sector of the matrices of differential op- 
erators behave the same way as c numbers, when the boundary data for 
integration is set to zero. This leads to the relation 



*i \ ^ / (ei)2(e2e3)-i -6263^ -636^^ \ ( '^^ \ 

^2 =2 -^1^3^ (62)^(6361)-^ -636^^ 

*3 / V -6162^ -626^^ (63)^(6162)"^ J \ V'S J 

The Gauss' law constraint must be solved in conjunction with the Hamilto- 
nian constraint, given by 



^Vari^ae) + Adet(*„e) = 0. (108) 

Due to the invariance of the quadratic part of (108) under rescaling, we will 
consider a special case. Define 
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^'a = mat''-''; ifa = n^e^^-^ 



(109) 



where nia and are spatially constant. ^'^ The the eigenvalue of the action 
of Ea and its inverse on the exponential in (109) by 



e„{e*^-l = (yl-i)^„ Ae*^-^ = e^e*^-; (110) 
where we have defined = {A~^)akj. Similarly, 

e-\e^-^^} = A« j dx^{e^-'^} = e^-He"^'''}, (HI) 

where we have defined = '^jAj/kj. Note that only the projection 
{A~^)l^ onto the wave vector ki enters into the Gauss' law constraint. 
The Gauss' law constraint then reduces to 



-6362 ^ 



It is convenient for calculational purposes to define Gae — 

(eb) ^, then the 

map from the anisotropy to the shear elements is given by 



*3 = ^ (-ei2¥'i - e2i(/?2 + e3ie32(/?3) ; 

*i = ^(ei2ei3¥'i - e23¥'2 - 632(^3); 
*2 = ^ {-ei3(pi + 623621(^2 - esiifs) . (112) 
The Hamiltonian constraint is given by 



^-3 ^ 




-6162 ^ 


-6265"^ 


^1 


-e?(e2e3)"^ 


-626^^ 


*2 J 






-eKesei) 





f ^1 \ 


II 









= -For* + Adet^ = 0. (113) 

The ingredients for (113) in terms of the CD J matrix elements are given, 
starting with the variance, by 

^^Note that this an example of a configuration where Af are spatially homogeneous, 
but are not. 
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where we have defined |ep = (ei)^ + (62)^ + (63)^. The determinant is given 

by 

(115) 

We will show the intermediate stages of the calculation of (112). First we 
have 

2^^1*2*3 + (pi<P2(P3 = ^{A + ipi(p2(P3) ; 

+ V32*i + V'3*i = \{A + 6^1<^2<^3) , (116) 

where A is given by 

A = (ei2ei3)^v'-J - (ei3)^(/92V?iV92 - {euf ^pi(pnp3 
+(e23e2i)^V'2 - (e23)v'2V'2V5i - (e23)^</'29'2V'l 
+ (e3ie32)^V'3 - (e32)^</'3</'3V'l - (e3l)^V'3V'3V'2- (117) 

The main observation is that the difference between the first and the second 
line of (115), which is det^', vanishes. By virtue of the choice of configu- 
ration the CDJ matrix is degenerate, and is insensitive to the cosmological 
constant. The Hamiltonian constraint (113) then reduces to the variance, 
which is the condition that (114) be zero. This is simply a quadratic equa- 
tion in (^3 with roots 









[(!)■ 













Substitution of (118) into (112) yields the shear CDJ matrix elements 

*i = ^( -— (1 - — - U + —)<P2T—V^T^]; 

2 V 62 V 63/ \ 63/ 63^ J 

^2 = l(-(l + -)m--(l- -)^2 T -V^) ■ (119) 

2 V V 63 / 61 V 63 / 63 J 
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Since the CDJ matrix is degenerate, then the 3-metric hij cannot be con- 
structed. Nevertheless, the solution for is well-defined and provided in 
terms of (pi and ip2, regarded as the two physical degrees of freedom per 
point. The information regarding the connection is encoded in the coeffi- 
cients Ca, and is coordinate-independent. 
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4.2 Passive difFeomorphism configurations 

The previous example is a special case of the following situation where the 
connection Af is not spatially homogeneous, but corresponds to the Jacobian 
matrix of a transformation between two coordinate systems. For each point 
P in the 3-dimensional spatial manifold S, there exists a homeomorphism 
from S to which expresses the coordinates of P in a local chart. The 
atlas of S consists of the set of all local charts of S. We will associate each 
chart to the Ashtekar configuration space variables as follows. Define an 
invertible transformation from a chart depicting local coordinates into a 
chart depicting local coordinates (f)^, where 

(1)^ = (l>\x,y,zy, = (t>^{x,y,z)- cff' = (tf'{x,y,z). (120) 
The transition functions define a Jacobian matrix J, given by 




(121) 



and for J invertible, the map x (/){x) constitutes a homeomorphism. Now 
choose the connection to be the Jacobian matrix 




(122) 



The choice (122) can be thought of as a triple of flat abclian gauge fields. 
However, the magnetic field is nonvanishing since it is a nonabelian gauge 
potential. The Ashtekar magnetic field is given by 

where the abelian part vanishes and we are left with the nonabelian part. 
The helicity density matrix is given by 




(124) 



Since Cae is isotropic, this configuration is definitely integrable. The vector 
fields are given by 
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which are exactly invertible since they are variations with respect to local 
coordinates. 

To show that the structure functions needn't be numerical constants, let 
us compute the algebra of vector fields 



[Va,Vb] = /if^V/. 



(126) 



Expanding this we obtain 







/dx-> ^ 


— 1 




)-, 

' dx} ' 




dxU 



d'^x^ 



dx' '6^x^ 



\ d 



/ \ / d(p^ \ r dx' 
From (127) we read off the structure functions as 



(127) 



hi.dx'^ A dx' = (1^) (1^) {dx\ d{dxydcl>^)} 



(128) 



which take on the interpretation of the Poisson bracket between the coordi- 
nate and its velocities with respect to the new coordinate system. 

The Gauss' law constraint for this class of configurations is given in 
matrix form by 









d 





d 



*3 ^ 




( 


a 


^1 


V 







*2 / 




\ 








d 
dp 






Note that the magnetic hclicity density contribution is zero, hence the in- 
tegration of this system should be relatively simple. To solve this it is 
expedient to map into the corresponding coordinate system 



<Pf{x,y,z)^ipf{(P\(t>\cP^), (129) 

since the coordinates t are globally well-defined within the domain of defi- 
nition of the passive diffeomorphism. The solution is then given by 







(s) 









d 



dd 







d 

w 



\ -1 / ^ 





A / 





d 








V 
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The cyclic operators are given by 



7r2 = 82 ^dsd^ ^^2^3 ^di = 1; 
TTs = 8^^8182^838^^82 = 1; 

TTl = 8^^828^^8182^83 = 1, (130) 

which are the identity since the vector fields and their inverses take on the 
clear interpretation of differentiation and antidifferentation with respect to 
the chosen coordinate system. The intermediate stages of inversion of the 
shear subpsace kinetic operator is given by 




92"' \ / 1 838^^828^' -838^' 

8^^ ^ -^182^ 1 8182^838^ 

5f ^ / V 828^^8182^ -8283^ 1 



82' 


\ 




81^82 - 


-8381 







1 f 'did^' 


1 


82^83 





8^') 


V d^'8i 




1 




O2' 




-82^838^ 




K 


-8^^8182' 


93' 


82' 


') 




93' 


-8^^828^' 


9i' 





Prom this we can construct the Gauss' law propagator 



/ _a-ia2 {83?{8l82)-^ \ 

^1 = ^ {9i)\8283)-^ -828^' -8382' 

\ -818^^ {d2?{d3dir^ -838^^ I 

and the general soluution to the Gauss' law constraint is given by 



^ae = ne^g = Tae[0:% T^e = ^^le + EL-Jfil]- (131) 

The Hamiltonian constraint can then be written as a condition on the re- 
duced space of diagonal CDJ matrix elements as 



Qf'3'<pf,ipg, + Qf'3''''ipfipg>iph' = 0, (132) 

^^Equation (131) and what follows can be taken as the general notation for all cases, 
not limited to the one being considered in this example. 
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where we have defined the following operators from the invariant tensors of 
the unconstrained space ^oe, 



Qf'9' = 1 (^«e^6/ _ S''fS''^)T£T^'f (133) 
from the variance, and 



Qf'9'h' ^ -e-'^-e^f^T^Tf^T^; (134) 

from the determinant. Equations (133) and (134) can be seen as embed- 
ding maps from the physical degrees of freedom, which are independent of 
coordinates, into the bigger space which takes into account the kinematic 
effects of the Gauss' law constraint. The inverse map constitutes the projec- 
tion to these physical degrees of freedom, of which two per point are freely 
specifiable. This can be realized when one solved (131) to yield 

</f3 = (/?3(¥'i,(^2;7)- (135) 

The notation in (135) indicates that each solution can be expressed in terms 
of two degrees of freedom per point ipi and ip2, in conjunction with the 
equivalence class of integral curves 7 which completely define the Gauss' 
law constraint. The spatial 3- metric can be reconstructed from the solution 
to the initial value constraints via the relation 

hij = r%B-')t{B-'rj{detB), (136) 
where we have defined the internal metric 77"^ by 

r?- = (det*)(*-i*-i)- = idetT[^]){T-^[0\T-'[0\r^ = r?-[<^i, <^2; 7](137) 
The metric is given by 



hij = ^""-Q^Q^ = hiM^ '^2;7]- (138) 

One could present the argument that all magnetic fields within the class 
of passive diffeomorphism configurations fall within the same equivalence 

class of the constant field = 5ia, since they are obtainable by coordi- 
nate transformation. However, the metric still has a complicated nonlinear 
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dependence on the particular within this class, owing to the nonlin- 
ear Hamiltonian constraint. Equation (132) can be solved formally using 
Feynman-diagramattic methods, which we will not display here. 

The previous example treated is a special case of the passive diffeomor- 
phism configuration. Corresponding to each general coordinate transforma- 
tion is a configuration solving the Gauss' law constraint as indicated. To 
complete the solution we incorporate the Hamiltonian constraint. 
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4.3 t'Hooft Polyakov monopole 



We will now consider configurations Af which do not correspond to a well- 
defined map between two charts, for which the Gauss' law constraint is still 
integrable. Let us examine the case of a nonabelian magnetic monopole 

^" = ^-i^' (139) 

of monopole strength g, where in this example and the remaining examples 
to follow, = x"^ + y'^ + The magnetic field is given by 

and the helicity density matrix Cae is given by 

Cae = ^''^^~'^\ earjX^x'x') = 0, (141) 

which vanishes due to antisymetry of eijk- Therefore we should expect this 
configuration to be integrable. 

The magnetic field can be decomposed into a product i?* = Rj(6jaB^), 
where 



i 5 1 / \ 

This suggests that oc (x")^ is the part of not constrained by coordi- 
nate system. The part intertwined with the coordinate system is i?*-. Note 
that dct(Rj) = 0, causing to be degenerate. Nevertheless, we can still 
solve the constraints. 

Note from (140) that Na oc X a, yielding a vector field 



v„ = X„ + iV„ = ^i^x«x^A. (142) 

Since Xa is integrable, and Na oc Xa, it follows that Na and therefore 
Vq is integrable. Indeed, defining Va = x^x^di, one sees that [fa,^^6] = 
with vanishing structure functions /i^^. The integral curves of = ^, 
correspond to a flow 



^^This can also be seen as an instanton of zero size. 
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(143) 



Cancelling off a factor of r ^ from (142) and using Cae = 0, the Gauss law 
constraint for g ^ 2 reduces to the form 

z y \ / d^i/dt \ / X \ / dcpi/dt \ 

2 x d^2/dt = - y dip2/dt . 
y X J \ d^s/dt J \0 z J \ dips/dt J 

We have defined a parameter t to denote location along the integral curves 
of the vector field x^di, which acts as the identity operator on coordinates 
^ = X* . Using the relation 

we can reduce the Gauss' law constraint to the form 



Eliminating the differential operator ^ — 1 from both sides, we can a 
braically solve for '^a = ^a['p] 

/ z y Y' / X \ / ipi\ 
^3 J \y X J \0 z J \ ip3 J 



/ 


yz 


— K 

z 








X 

z 


vL 

zx 


i] 




V 


x_ 


— U. 




V <^3 


y 


X 


xy / 



Expanding this out we have 
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In conjunction with solving the Gauss' law constraint, one must also solve 
the Hamiltonian constraint in order to obtain a full solution to the initial 
value constraints. The Hamiltonian constraint is given by 

+K[lPHP2(P2, + 2*1*2*3 - <^i*? - <^2*i - <^3*i) = 0. (146) 

Since the relation between * and (p is algebraic, then the Hamiltonian con- 
straint can be solved exactly in closed form. We will do this, by substituting 
(145) into (146). The intermediate stages of the calculation are as follows. 
For the shear contribution to the variance we have 



+2 



2 y2 



fz^ — x^ — y'^\ fx 

\ i2 



2 2 2 2 

X 2 (y r \ 2 
z'^\ (y^ — z^ — x^ 

— J^2<^3+ ( -2 



¥'3</Cl 



(147) 



The shear contribution to the determinant is given, after a long algebraic 
calculation, by 



*1*2*3 = \ 



X \ 2 (y \ 2 



^j^2^3-(^jv'2V^3-(-2jW-(-2 



.z^' vz^/ \z^ ' vx^/ ^y ' ^y 

The other contribution to the determinant besides (pi<f2V3 is given by 



(148) 



\y2^2 J 



A 



,2 , 



(149) 



We have shown the intermediate stages of the calculation to point out a 
remarkable effect which happens for this configuration. Upon insertion of 
the results of (147), (148) and (149) into (146), the coefficient of A vanishes 
due to a precise cancellation and the Hamiltonian constraint reduces from 
a cubic to the following quadratic equation 
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+ V'2</'3 + - {'^l? - {"^2? - (^3)^ = 0. (150) 



The result is that the monopole configuration is insensitive to the presence 
of a cosmological constant. Putting (145) into the right hand side of (150) 
we obtain the following quadratic equation 



^2^3 ^3^1 



0. 



(151) 



x^ 

Equation (151) can then be solved for ip^ = (p3{(pi,(p2) and the remaining 
CD J matrix elements determined. The result is given by 




■^2 = - (^) VI ± (^) Vv^; *3 = ±x/^m, (152) 

whence *ae has been reduced to two physical D.O.F. While 5* and ^'ae 
are degenerate as is the 3-metric hij, the CD J matrix still has a solution as 
in (152). One may attempt to determine the algebraic classification of the 
spacetime that this describes as in Paper VIII. The characteristic equation 
for the CDJ matrix is given by 



- (tr^')A2 + ^(Far*)A - det* = 0. (153) 

The invariants of "^ae can be computed from (152), as well as what we have 
determined on the basis of the Hamiltonian constraint 



Var^ = deW = 0; tr^- = ipi(l + ^^+ (l + ^^± -^^flpwi- (154) 

This reduces (153) to A^(A — tr^*) = 0, with a double root A = and a single 
root A = tr*. A degenerate CDJ matrix with two equal eigenvalues implies 
a spacetime of Petrov type 111. 
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4.4 Another extreme example 

Let us now examine a potential with a different structure, 

At=g^. (155) 
Tfie magnetic field for (155) has an abelian part 

— r) 

and a vanishing nonabelian part 

(BlUna^eHan = ^^^^^ V'^''^ (^) (^) = 0. (157) 

This is an example where Xa 7^ and Na = 0, which by the previous 
analyses is integrable.^^ The helicity density matrix Cae also vanishes, since 



Cae = AiBi = ^e'^'^x^x'x^ = 0. (158) 

Note that Af and have switched roles in relation to the monopole example 
treated earlier. The Gauss' law vector fields with rescaled versions are given 
respectively by 



These vector fields satisfy the commutator bracket [va,^^] = eabcVc for angu- 
lar momentum, with numerically constant structure functions. The Gauss' 
Law constraint reduces to 




Vi 


















II 










^^3 / 




Kf3 J 



One may solve for ^'y = 5'/ [(^ as 

^ = {l + Tr)-^v'^Mip = J^ipg, (160) 



'^"Thc magnetic field is degenerate for this configuration, which malces the 3- metric hij 
also degenerate. 
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where we have defined 



as well as 




V2 
i;3 I ; TT 
vi 



-V3V1 ^ 

M = I —V1V2'' 1 ViV2^VsVi^ 




The cyclic operators ir^i are given by 



7r2 = V2 ^V^Vi ^^2^3 ''"'^1; TTa = ^ViV2 ^^3^1 ^V2; TTi = ^V2V^ ''"^1^2 """f^S^l^l) 

With all the constituents of the inversion defined, we must define the action 
of Va and in local coordinates. Since Va generates rotations in the plane 
orthogonal to the coordinate Xa, we coordinatize the integral curves using 
angles ^q, for a = 1, 2, 3, such that 

= ^ ^ = J d4>a. (162) 

Note that x"" is the spectator variable in (162). Along with the identifica- 
tions (162) comes a prescription for changing into the appropriate coordinate 
system to carry out the inversion. 

When acting with V'^ and v^^, one performs the transformation (x, y, z) — >■ 
(r, ^3,2;) and its inverse, which leaves the z direction invariant. 

X = A/r^ — z'^cosct)^; y = \J r'^ — z^sinc/is; z = z\ 

r = \/a;2 + j/^ + z^; ^3 = tan~^ ^— j ; z = z. (163) 

When acting with vi and v^^, one perform the transformation {x,y,z) 
(r, z) and its inverse, which leaves the x direction invariant. 

y = \J — x^cos(?!)i; z = \J — x^sin^i; x = x; 

r = \/x'^ + + z^; 4>i = tan~^ (-] ; x = x. (164) 

\y/ 

When acting with V2 and f^^, one perform the transformation {x,y,z) — >■ 
(r, (f)2, z) and its inverse, which leaves the y direction invariant. 
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r = ^/x"^ + + z^; (p2 = tan~^ (^^y, y = y. (165) 

The end result is that the shear elements are images of anisotropy ele- 
ments (ff under coordinate transformations. Hence the latter play the role 
of the (three) physical degrees of freedom, with the former unphysical. In 
this example, both Af and are degenerate, whic means that the spatial 
3-metric hij will be degenerate in this case. 
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4.5 Instanton of nonzero size for g — 2 



To provide a nondegenerate integrable example where Cae is nonvanishing, 
let us consider the case of an instanton of size p and strength g. We will 

approach the problem in stages, first considering some simple cases and then 
moving on to the general case, to illustrate some techniques in integrability. 
Starting from a connection of 



= ^rr^^aijx^^ (166) 

the magnetic field and helicity density matrix are given by 



= {^.2^ p2y {-'^9P^5ai + gig - 2)x*Xa); Cae = -2j-^^^^^eaejX%167) 

For this configuration the Gauss' law constraint is given by 



2 r/ vy v,\ / 3y 3z \ 1 / *3 \ 



+ 



(r^ + p2)2 



Vx 

Vy 



+ 





where we have defined 



j,2 _|_ p2 



2x —X —X 

-y 2y -y 
—z —z 2z 



"Pi 

<P2 



= 0, 



d / g — 2\ d d d d d 
and similarly for Vx and Vz-^'^ The differential operators satisfy the algebra 



[dt^d^i] = -d^i Vi (169) 

We will first treat the g = 2 case, since the contribution to (168) proportional 
to dt vanishes which simplifies the inversion. We will retain the symbol g 
in what follows, with the understanding that g = 2. Note that the helicity 
density matrix Cae can be combined with the spatial gradient terms using 
an integrating factor 

^'^ Recall from previously that dt is integrable. 
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This is possible since Cn = C22 = C'as = 0, which corresponds to one of 
the off-diagonal helicity configurations, whence (25) reduces to just the first 
term. For g = 2, the Gauss' law constraint reduces to 

dy dx 



da-. 



-X —X 



ir' + pY'i -y dy -y \{r^ + p^r 



-z —z d. 



z 




Then the solution is given by 

= H^f. (171) 
where the Gauss' law propagator is explicitly given by 



dy <d dz\ I dx 



H = -{r'+p')-''" \ {r'+py -y dy -y {r'+p 




dy dx / \ ~^ ~^ ^; 

The operator ordering must be preserved, but it is straightforward to carry 
out the operations, which have been reduced to computing the inverse 

dy d,\ 

dx = {l + Tr)-^v-^M. 

dy dx J 

Using the results from the previous section on inversion, one defines 



dy 



■'z 

00a 



TT 



X 




as well as 

1 dzd~^dyd^^ —dzd^^ 

M = I -dxd-^ 1 dxdy^d.d-^ 

dyd-^d^d-^ -dyd-^ 1 

The cyclic operators n^i are given by 
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-Kz = d^^dxd:^^dzd~^dy; 



TT. 



X - d-^dyd-^dxdy^d,. (172) 

If one restricts oneself to integrations wlicre the integration constants are 
set to zero, then the operators dx and d^^ are inverses of each other. All 
operations in the cyclic operators cancel, yielding tTj; = tTj^ = tt^ = 1 which is 
the identity operator on functions, implying that (1 +tt)~^ = ^. Analogous 
cancellations occur in M, and we are led to the following result 



1 d^ ^dy -d,d^ 

Vae = \ dx d, ] II -dxd-^ 1_ d-^d, 

d^ ^dx -dyd^ 1 1 

5-1 5-1 -5-15,5-1 

-5-15,5-1 5ji 5-1 
5ji -5-15,5,-1 5-1 

One can verify by the multiplication of the operator matrices, that v^^v^f = 
5a f- The end result for 5 = 2 is that one can define the Gauss' law propagator 
Tg as an operator-valued matrix, such that 



dy 





5. 


dx 


5. 








dy 


dx 






( 




-w 





5-1 

Uy 








5-1 










^f=ff^g, (173) 



where 




Tl = -{r' + P'r''/' -dz'dxdy^ 5-1 

V 5-1 -5-15,5-1 

-y dy -y ]{r' + p'y, 

-z -z dz J 

The operator ordering must be strictly maintained in this inversion. Accord- 
ing to (173), the shear elements ^ f are a map of the anisotropy elements ipf 
under the action of Jj. There is a reduction in by three D.O.F., conse- 
quently we for the purposes of this paper we treat as the physical degrees 
of freedom. The action is well-defined, as one can verify by multiplication 
of the matrix elements. The result is that the metric for an instanton on 
nonzero size should be nondegenerate. 
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We have shown that the Gauss' law constraint has a well-defined solution 
for g = 2 for the instanton, but we would like to extend it to more general 
values of g. One possibility is to modify (173) to 



*/ = T>V,-p^*„ (174) 

where we have defined 

\ y X J 

which parametrizes the deviation from g = 2. The operator dt has a simple 
action on polynomials, and for ^jv(r) € Hn, the space of homogeneous 
polynomials, given by 

An{^= Mmnx'y'^Z^, (175) 

we have that 



^^A^{x) = N A^{x); d^^A^{x) = ^AN{x). (176) 

The general solution to the Gauss' law constraint for this configuration is 
given by 



= (177) 
where we have defined the Gauss' law propagator J| by the infinite series 

oo 

J] = {sj + =sj+ j:i-^rpfpf, . . . PtA.- (178) 

n=l 

A necessary condition for convergence of (178) is that — 2| < 2p^. 

We can get a better idea of the action of the differential operators by 
explicitly evaluating their action on a set of basis functions, first considering 
the effect without the magnetic helicity density. Take Vx for example, where 
we wish to solve the differential equation 



d d \ 

d-x-'^dir='^- 



(179) 



47 



Equation (179) yields 



1 — k J dxx-^ 



Q, 



(180) 



which upon choosing Q to be any single term of the form (176) expands into 
the infinite series 



^ = J2[k J dxx-) {x'y^z^). 



(181) 



The zeroth order action is the identity. The first order term is given by 



k j dxx^^{x^y"'z'') = k{l + m + n) j cZa:(x'+ V^") = fcx^ (^^^i^) {x^ y"^ z^%\^2) 
The repeated action leads to the following series 



+k\ 



,,l + m + n\ 9 , 9 + m + n\ // + m + n + 2 
l + k{ — ; — - — ]x'^ + k'^ 



1+2 J '"V 1+2 J\ 1+4 
l + m + n\/l + m + n + 2\/l + m + n + 4: 



1+2 J\ 1+4 
This is a hypergeometric series 



1 + 6 



Kl + m + n + 2\ f i + m + n + 4:\ 
1 + 4 )[ — TTe — > 



.(183) 



^ = X 



[0, ; -; kx^j , 



and the ratio test yields a radius of convergence 



(184) 



<x < 



k 



ff-2 
2fP ■ 



(185) 



Performing the analogue of (180) for the y and z directions enables one to 
construct the Gauss' law propagators Jj and the 'deformed' basis elements 
Tae- Provided that the physical degrees of freedom are restricted to the 
basis functions (176), then one may solve the Hamiltonian constraint (132) 
with respect to this basis. Since our main purpose in this paper is to demon- 
strate the integr ability of the Gauss' law constraint, we will not carry out 
this exercise here. 
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4.6 Instanton of nonzero size for 5' 7^ 2 



We will now treat the general case of the previous example, using the general 
technique of integral curves. First note that the magnetic field can be written 
in the form 

^ / ^ + (mx)^ m?xy 

^ \ m zx m zy 

where we have defined the mass scale 



1 

m = - 

P 

We can directly read off the vector fields as 



m xz 
m^yz 
1 + {mzY 




Wx = (1 + (mx) )— + {mx)[y— + z-^J + 
Wy = {l + [myf) ^ + (my) (z^ + x^ + 



dy dzJ'p'^ + r^^ 
d / s / 5 d \ 3g 
dy \ dz dxJ p'^ + r"^^ 

(1 + imzf) I + imz) {x^ + 4) + TTT^- ^'''^ 



Since is conformally related to a magnetic field of the correct form, 
this means that the system is integrable in terms of coordinates. We will 
illustrate for the x direction with the remaining cases to follow by cyclic 
permuatation. We will take the integral curves as 7 = {u,v,w). From (187) 
we read off the following relations 

^ = 1 + (mx)^; — = (rnx)y; ^ = (mx)z. (188) 

au au an 

The first equation of (188) integrates to 



u — uq= I -. — )■ X = — tan [tan ^ {mxo) + miu — uqVi , (189) 

J 1 + {mxY m '- ^ ' 

which defines the projection of the integral curve u onto the x direction. To 
find the projections onto the y and z directions we must substitute (189) 
into the second and third equations of (188). For the y direction this yields 



— = {mx)du = tan [tan ^(mxo) + m{u — uq)] du, (190) 
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which restricts the x dependence to the curve u, integrating to 



y= ^2 ; fO (.191) 

cos [tan ^ {tuxq) + m{u — UQ)y cos [tan ^ (rnxo) + m{u — uq)] ' 

We have included the analogous steps for the y and z directions. 
Let us define the following functions 



ri{u;xo) = tan [tan ^(mxo) + m{u — uq); 

({u; Xo) = cos [tan~''"(ma;o) + m{u — uq)] . (192) 

Then the integral curves 7 are well-defined/^ and given by 

xiu) = ^-^; yiu) = ^^y y(-) = ^(^ (193) 
for the u curve, 

C{v;yo) m C{v;yo) 

for the V curve, and 

x[w) = — r; yiw) = — -; z{w) = (195) 

C(w;zo) C{'w;zo) m 

for the w curve. Having constructed the integral curves of the vector fields 
Va, we must now take into account the helicity density contributions. First, 
note that (187) simplifies to 



^ du p^ + r'^{u)^ ^ dv p^ + r'^{v)^ ^ dw p^ + r'^{w)' 

We can now treat incorporate the helicity density contributions as integrat- 
ing factors, restricted to the appropriate integral curves. Make the following 
definitions for integrating factors 

, r rr 2 , ( r]{u'-xo) ^^ ^ ( yo ^\ f zp ^2^-1^ , 

^„ = exp[,y^jp+(— ^) +(^^) +(^^)J du; 

0. = exp[,yj, +(— ^) +(^^) +(^^)J dv; 

= expL r + (!^)^ + (-^y + (-^)r^Mi97) 

I- 7«,o \ m J \Ciw';zo)y \C{w';zo)y J 

^^Tho parts where r/ and blow up correspond to infinite values of x, y and z, therefore 
these integral curves should cover E witfiin a finite range of the parameters of the curves. 
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While the results of the integrations are unwieldy to display in closed form, 
the point is that the integrating factors are path independent and well- 
defined for g <2. 

We can now write the inverse Gauss' law operator on the shear subspace, 
rigorously, as 



-1 



w 



•yOdyO(f)^^ (t>W ° dyj O (f)^ 

1 _ I J. ^ « ^ ^-1 ^ ^ « ^ n 1 _ n_L^^-l»..-l 



0«oa„o0^i o dyj o (f>^^ I = (I+tt) 

(pvodyo (p-^ 4>u °duO <t)-^ 



where we have defined 

/ 1 (f)uO o (f)uv odyo -(/)^ odwo 

M = -(pu oduO (puv o 3^7^ o 1 (j)^o 

\ (t>w° ° (pwu oduO -0^ odyO o o 0-1 

with (j)uy = 0-1 o 0^, and 

/ 0^ o o 0-1 \ 

w = { (pyjodyjo 0-1 j ; vr = 

V 0« o 5„ o 0-1 y 

as well as the cyclic operators 

vr2 = TTa = TTi = 1. (198) 



7T2 






















TTl 
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4.7 Another nondegenerate example 

For a potential of the form Af = Sai{g/r) for r > 0, the magnetic field is 
given by 



K = Sa^{^)\e,aJX^[^), (199) 

and the helicity density matrix by 



Cae = AfBi = 6ae[^) + EaejX^ {^) , (200) 

which is an off-diagonal helicity configuration. The abelian part of yields 
a vector field 



^« = ^a^^ = {^)^a^JX^Q^ = {^)ja, (201) 

where Ja are the orbital angular momentum operators, which satisfy the Lie 
algebra [Ja, Jb] = i^abcJc- The integral curves of Jj correspond to a flow 



(t)0^ {x'^,x^ ,x^) = [x\ x^cosdi - x'^sinOi, x^cosOi + xhmBi) , (202) 

which leaves invariants with {x\x^,x^) = Perm{x,y,z). Note that Jj 
annihilates functions of r as well as functions of the single variable x*. The 
nonabelian part of Va is given by 



Nn. 



2 d 



dx"- 



The Gauss' law constraint is given by 



(203) 



Vy 





Vz 


Vx 


Vz 








Vy 


Vx 



+ 





+ 



Vx 











Vy 











Vz 



2x 


— X 






-y 




--;)] 


(i) 


—z 


—z 


2z J. 





0, 



where we have defined Vx = Xx + Nx, with the Abelian and nonabelian 
contributions are given, up to a factor of (g/r)'^ by 



52 



(^)j^; N^ = d^ (204) 



*1 \ 






X 


X 






I ? 


Vy-^ 


_ y 


^-3 J 






z 


Vz-^ 







H 









and similarly for Vy and t;^. 

Since Cae is an off-diagonal helicity configuration, the Gauss' law con- 
straint can be written 



Vy V, 

14 v; 1 I I + I v,-^ I I I = 



where we have defined 



yx = ^ + — + — 205 

ox rg' 

and likewise for Vy and V^. The solution is formally given by 

= j/<^/, (206) 
where we have defined the Gauss' law propagator 



Vy Q v. ' ^ 



H = -\ V, 



Vy 



T/ ^ 


a; 


X 


— 4- 




'y 
--pi 


z 


z 


Vz-fi 



for the chosen configuration. Using the results from the previous section on 
inversion, one can explicitly write the inverse 

fvy v,y' 

^ae = = (1 + -kY^-^M, 

V Vy V^ j 

where the constituent parts are given by 

1 v.v-^VyV-^ -v,v-' 

M = I -V^V-^ 1 V^V-^V,V-^ 



and 



VyV-W^Vy-^ -VyV-^ 



Vy ^ ^ 








-Ky 











T^z 











T^x 
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as well as the cyclic operators 



% = v-^v,v-^VyV-^v,- 

T^x = V-WyV-W^Vy-^V,. (207) 

To give this all meaning, one must understand the action of V~^. Let us 
demonstrate with the x direction for illustrative purposes. To find the action 
of is the same at attempting to solve the equation 



^ ( A + ^ + ^(^, = (208) 

\ox rg/ 

for ^ in terms of Q, given ^ as the image of Q under some map. Note 
that Jx = ydz — zdy has a trivial action on the x and the r dependence of 
operations, since 

Jx{x} = Jy{y} = Mz} = J^{r} = Jy{r} = Mr} = 0. (209) 

We may exploit this observation to construct an integrating factor operator 
corresponding to each direction 

(l)^ = r^{2{r + x))a'^''; (l)y = r^{2{r + y))a'^^ ■ (P, = r^{2{r + z))o'^% (210) 

where the coordinates are rescaled to make them dimensionless. This trans- 
forms (208) into the form 



0x ^dxi-^xft) =Q — >^ = (f>x^ J dx(f)xQ(,x, y, z) 



(211) 



which integrates the x dependence while holding y and z constant, and 
similarly for the other directions. It is convenient to think of this as a 
composition 



= </>-l o a, o = <^-l o o (212) 

taking care to preserve the operator ordering. Then (j)x is nothing other than 
a homeomorphism from S onto a particular chart adapted to operations in 
the X direction. Define transition function operators Uy.ixj = 4>x^'P~^l- Then 
the cyclic operators (207) are given by 
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T^y = (Py^ ° O Uyz od^o Uzx ° O Uxy o Oy O Uy:, o d^^ o U-^x ° dx o (f)x 

Tfz = o o Uzx odxO Uxy o dy^ o Uyz odzO Uzx o o o dy o 

= (kx^ ° dx^ ° Uxy odyO Uyz O O Uzx °dxO Uxy O O Uyz O 8^0(1)^ 

(213) 

If not for the the transition operators Uxy, the cychc operators tTx, T^y and 
TTz would be the identity. Indeed in the strong couphng hmit — )■ oo, we 
have that 



hmg^oo^x* = \vcQ.g^ooUxixj = 1- (214) 

So the transition functions Uxixi = ^x^x^id) parametrized by g.^^ The 
matrix M is can be written as 

/ 1 VyOTTyOVx-' -V^Vx-' 

M = -VxV-^ 1 K O TT, O 

or in component form taking into account the composition operators 



Mae = ^ae " 5ay5ex4'x ^ ° dx o Uxy o dy ^ o (j)y 
+5ax5ey(ky^ O dy O (py O TTy O O O (px + . . . , (215) 

where the dots cycle through the cylic permutations of indices {x,y,z). 
Hence we have provided a prescription for inverting the Gauss' law constraint 
for this configuration, which leaves remaining the Hamtilonian constraint to 
fix the cpf. Once this is completed, then the 3- metric can be constructed 
from the CD J matrix and the magnetic field from (199), as in 



h,j = {detm)&ffVgVg,)-\B-')^{B-'y., (216) 
where we have defined 



- ae — ^ae ' ^ae'^ f 



+ eUI (217) 



^^In the large g limit the computation of j/ proceeds straightforwardly as in the previous 
subsection, so we do not repeat it here. One could envision performing an asympotic 
expantion in inverse powers of g to find the Gauss' law propagator for large g < oo. 
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4.8 Diagonal connection 



Finally, we will consider a diagonal configuration with three degrees of free- 
dom per point, for example one of the quantizable instanton representation 
configurations identified in Paper XIII. The connection is given by 

I A\ Q 
A" = 
\ 

This yields a magnetic field and helicity density matrix, respectively, of 



(^2^3 ~dzA^ ^2/^3 \ / ^1^2^3 ~-^dz^ ^l^J/^3 \ 

dzA\ ^3^1 —dx^i j ; Cae = I -^d^A^ A^A^A\ —A^dxAi^ 1 

We have chosen x, y and z for coordinates, without loss of generality since 
these are merely labels. The vector fields Vq are given by 



^ ^dx \ dy J dz \ dz / dy' 

V2 = A'A'-+(^)--(^)-; 

^ ^ dy \ dz ) dx \ dx ) dz^ 



V -4142 5 , d (^Al^ d 

V3 - A,A,- + - - [-^) -. (218) 

We see the familiar decomposition of Va into a part containing a Hamilto- 
nian vector field Xa and a part containing a nonabelian contribution Na- 

In the analysis of the previous section we have treated Xa and Na as sepa- 
rate, delineating the condition for each part individually. But we will show 
that the present example is an example where the two vector fields can be 
combined into one. 

For Vl the Hamiltonian function is A\, which induces a flow in the two 
dimensional hyperplane dual to the x direction according to the Hamiltonian 
equations 



du dz ' du dy 

Since A\ is assumed to be a smooth function of x, y and z, then there exists 
a Hamiltonian flow generated along a one dimensional curve 7a; (y, z) 



y = y{yo,ZQ;u) ; z = z{yo,zo;u) 



(220) 
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in the {y, z) plane for each fixed x, where u parametrizes location along the 
curve. More precisely, (220) is the projection of the congruence 7 onto the 
{y, z) plane. To obtain the three dimensional picture, let us 'unfreeze' the x 
dependence and set 



^=AlAl = f{yo,zo;x,u). (221) 

We have expressed A2 and in a chart adapted to n, which entails sub- 
stitution of (220) into the y and z dependence of these functions. We must 
invert (221) along 7, which entails an antidifferentiation with respect to the 
parameter u. If 7 forms a congruence, then this should be possible since 
there are no sclf-intcrscctions of the integral curves of by supposition.^*^ 
Equation (221) is an ODE, which we should expect to have a unique solu- 
tion. A way to see this is using Pickard's theorem in the theory of exterior 
differential systems. Namely, given the one form 0^ = dx — f{x, u)du G B? 
with coordinates {x,u), then there exists a one-to-one correspondence be- 
tween solutions to (221) and curves '~f : R ^ B? such that 7(0^) = and 
^{du) / 0. The result is that restricted to the curve 7, the vector field V2 is 
integrable. Another way to see this is that 



de^ = -d^cfdx Adu — ^ 6>^ A de^ = 0, (222) 

namely that the topological torsion of 0^ vanishes when restricted to 7^. 
Applying a similar analysis to V2 we have 



dv dx ' dv dz 

Since A2 is assumed to be a smooth function of x, y and z, then there exists 
a Hamiltonian flow generated along a one dimensional curve 'yy{z,x) 



y{zo,xo;u) ; x = x{zo,xo;v) 



(224) 



in the {z,x) plane for each fixed y, where v parametrizes location along 
the curve. Hence, (223) is the projection of the congruence 7 onto the 
{z,x) plane. To obtain the three dimensional picture, we 'unfreeze' the y 
dependence and set 



^ = AlA\=g{zo,xo;y,v). (225) 



^"This should be the case if the integral curves are topologically equivalent to the Carte- 
sian axes, namely derivable from them by continuous deformation. 
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Equation (225) implies the vanishing 7^,(0^) = of the one form 6^ restricted 
to the integral curve 7^,, where 6^ = dy — g{y,v)dv. Likewise, the topological 
torsion vanishes along this curve 

de'^ = -dygdy Adv — > O'^ A dO"^ = 0. (226) 
Applying a similar analysis to V3, we have 



dx _ dAl_ dy _ _dAl 
dw dy ' dw dx 

Since ^3 is assumed to be a smooth function of x, y and z, then there exists 
a Hamiltonian flow generated along a one dimensional curve ^z{x,y) 



X 



x{xo,yo;w) ; y = y{xo,yo;w) 



(228) 



in the (x, y) plane for each fixed z, where w parametrizes location along 
the curve. Hence, (228) is the projection of the congruence 7 onto the 
{x,y) plane. To obtain the three dimensional picture, we 'unfreeze' the z 
dependence and set 



— = A\Al = h{xo,yo;z,w). (229) 

Equation (229) implies the vanishing jwi^^) = of the one form 6^ re- 
stricted to the integral curve 7^, where 6^ = dz — g{y,v)dw. Likewise, the 
topological torsion vanishes along this curve 



de^ = -dMz Adw — >e^A de^ = o. (230) 

We have seen, as a consequence of the Hamiltonian flow induced by the 
abelian part of the vector fields, that the vector fields are actually invertible 
including the nonabelian part, when one restricts the inversion to the integral 
curves 7 = {'yu,lv,lw)- Hence assuming that 7 is topologically equivalent 
to the Cartesian cooridnate system, then we have 

^r-^{^{u,v,w)} 
v^^{-^{u,v,w)} 



du^{u, V, w) 
J dv^{u,v,w) 
■- J dw'^iu, V, w) 



v.w const. 



w,u const. 



u.v const. 



(231) 
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for a smooth function Equation (231) is an abuse of notation, since 
^ may not initially be expressed in the chart {u,v,w). More precisely we 
must perform the composition of each opration with the conversion into the 
appropriate chart as in 



We will take these (j) operations to be implied, and revert to the former abuse 
of notation for convenience. 

Having put in place the aformentioned structures, we can now solve the 
Gauss' law constraint for the quantizable configurations. This is given by 
matrix form by 

^+Ci3(v) 

i + C23(n) iy + C2iH 

~-^2^z-^l ~^l^z^2 d^JJ + ^2(^1^2 

where we have defined 





ci3 = A\dyAl + 2AldyAl; 0^2 = A\dyAl + 2AldyA\; 
C23 = AldyAl + 2Ald,Al, C21 = Ald,A\ + 2A\d,Al; 
C31 = AldyA\ + 2A\dyAl; C32 = Ap^Al + 2AldxAl (233) 

adapted to the appropriate parameters for the given integral curves. The 
main point that this is an off-diagonal helicity configuration, which reduces 
the inversion of the Gauss' law constraint to the inversion of operators of 
the form 



du du 

W2 = ^+c{v) = e-^<^^'^^^e^<^)'^^- 
dv dv 

W3 = — + c(w) = g-Z^W^'"' Ag/cW'^"'. (234) 
dw dw 

The point is that restricted to the corresponding integral curves, one can 
write down integrating factors in the appropriate chart. The inversion of 
the Gauss' constraint then takes place exactly as for the off-diagonal helicity 
configurations. 
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5 Inversion of the vector fields in the general case 



We have seen that the Gauss' law constraint is exactly integrable for a large 
class of configurations for which Xa and are directly related, which es- 
tablishes a well-defined map (p — )■ ^[(p\- But we also wold like to establish 
the map for the general case of configurations A".^^ The sohition to the 
Gauss' law constraint requires a choice of 5* G C°°(S) combined with a 
well-defined prescription for carrying out the inversion. Whatever the pre- 
scription chosen, it is necessary to evaluate the inverses vj"^, ^ and v^^ 
as part of the inversion process. Assuming that this is well-defined then 
one ideally would like to integrate along the Va direction using holonomic 
coordinates, which requires the projection of onto the axes of a specified 
holonomic coordinate system. 

Let us say that we wish to invert the vector field Vq. First start with 
the differential equation 



^fa'^{x\x^,x^) = Q{x\x\x^), (235) 

where ^[x\x^ ,x^) G Dom\^a} and x-', x*^) G Ran{va} and i ^ j ^ 

A:^^ and the vector field is given by^^ 



Since Va might not necessarily be holonomic, we must specify a holonomic 
coordinate and then project Va onto the x^ direction. This choice cor- 
responds to the specification of a particular element = Sj^'* which will 
be divided out to normalize so that the integration can be carried out in 
this direction. Define the ratios 




(237) 



where ^ Q?^ Hence, (235) reduces to 

^^The purpose of the Gauss' law constraint is to establish the map from the physical 
degrees of freedom (p, irrespective of the wcU-definedness of the shear elements. 

^^The following argument can be generalized to all possibilities by cyclic permutation 
of indices. 

^*We will often abstain from using the Einstein summation convention, since we would 
like to explicitly depict the relation of the coordinate being manipulated to the spectator 
coordinates. Hence the notation (a;',x^ , x*^) will correspond to a given permutation of 
the Cartesian coordinates (a;, y, «) which are holonomic. The indices {i,j,k) each label a 
specific coordinate, as opposed to a range of coordinates. 

■^''Note, since it might not be possible to satisfy this condition globally in S, one is free 
to choose a different -B^ 7^ and then redo the procedure. 
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Next, act with the inverse differential operator corresponding to x'^. Hence 
(dx)^^ = [ dx and similarly for y and for z. By integrating both sides of 
(238) in the x'^ direction, we obtain 




(239) 



where {x'' , x^ , ^q) G Ker{d/dx^} is the boundary data, which must be 
specified on the x*^ = hyperplane. The solution then is given by 



+i5:2(x<,.^;.',5')£'<'||^, (240) 

where we have defined a parallel propagator with respect to the x^ direction, 
given by 

k 

K^{x\x^;x',0 = P{exp[- d^'' [fitk{x\x^ ,0^^ + ^%{^\^' ,0^)]) 



(241) 



Here in (241), P corresponds to path ordering along a path projected onto 

the x^ direction, the direction of integration in the inversion of the vector 
field Vq. For example when integrating in the z direction, the action of 
is given by 



^{x,y,z) = w^^Q{x,y,z) 
= exp (x, ^ + (x, 2/, ^) ) 

[*^-'°)^r^^'i^]- ''''' 

We have omitted the path-ordering symbol in (235), which should be un- 
derstood to act on both the exponential as well as the integral in the third 
line. 
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5.1 Gauss' law dual translation operators 



It will be convenient, for a certain class of configurations, to interpret the 
inverted Gauss' law operators as translation operators restricted to a two 
dimensional spatial hyperplane dual to a particular direction of path-ordered 
integration. The Baker-Campbell-Hausdorff formula for two operators A 
and B is given by [ ] 



exp(^)exp(5) = exp/x(A, B), (243) 
where iJ.{A, B) is given by 

H{A, B) = A + B + \[A,B] + ^[A, [A, B]] + 1 [5, [S, A]] 

+^ [S, [5, [B,A\]]+^^ [A, [A, [A, B]]\+... (244) 

The left hand side of (243) presents two disentangled operators A and B, 
which the right hand side combines into a single operator C = l-t{A, B). The 
decomposition into A and B is in general not unique, and may in general take 
on an intricate form due to noncommutativity of the operators. However 
for a given operator A, the operator B is uniquely determined once C is 
specified. 

To illustrate, let us express in the form (243) 



K,«(x^x^■;x^^') = exp - / dxHlk{x\x^.x^) 



exp 



dx^b%ix\x^ ,x'' 



d 
dx^ ■ 



where 6"^. and 6"^ must be determined. The quantities and ^5"^ in (241) 
are determined by S*, which fixes 6"^ and bjj^?^ For the operator ordering 

chosen in (245) , the action of is that of a sequence of translations in S 
by predetermined magnetic field-dependent amounts. Hence, 



(245) 



K^Q{x\x^,x'') = q(x'- j dx%1k{x\x\x^),x^ - j dx%%{x\x^ ,x^),x^^ 

(246) 

with cyclic permutation of indices, where we use the indefinite integration 
operator to signify integration in the x^ direction. Note that the x^ argu- 
ment of Q is unaffected by this integration, since that particular x^ depen- 
dence- is annihilated In* the d/dx^ and d/dx^ derivatives. The effect is to 

^"Onc imposes the condition that the left hand side of (245) decomposes into the right 
hand side with the operators djdx" and d/dx^ disentangled, and then use this condition 
to find the coefficients and Uj^- Alternatively, one may first choose 6"^. and 6"^, using 
this to define the ^f^. and ^"j, , which in turn fix the ratios of the Ashtekar magnetic field 
components. 
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translate the source Q according to 6"^ and 6"^., which are related to 
The are specified by choosing a well-defined configuration for i?*, then 
the procedure must be repeated for each B\. Note that the dual translations 
in the hyperplane commute with integration normal to the plane, 

hence (246). 

The procedure to invert the the Gauss' law covariant differential opera- 
tors is then as follows. First the practitioner must choose a configuration, 
namely a specified function = -B*(r) G C°°(S), and then evaluate the 
functional form of '^ae = ^ae[^6(^] for that particular configuration, for 
each configuration.^^ One chooses the six functions b^j^{r) and h":^.{r), which 
define a flow lines along which the source Qa becomes translated. One may 
then construct the Ashtekar magnetic field by first computing (3 from h using 
(243), (244),^^ hence the parametrization 

( Bl Q \ / HUb] PlM 1 \ 

5;= PUb] 1 

V 5f / V ^Ub] PUh] 1 / 

corresponding to an integration in the = z direction.^^ The Ashtekar 
magnetic field _B* defines three directions in S, one direction for each a. 
Hence, provided the Gauss' law constraint is satisfied for 9 = 9[\,Ba\, then 
the resulting Ashtekar curvature F^^ = Cij^B"^ constitutes a general solution 
to the Einstein equations by construction.^^. The physical interpretation is 
that Bl defines the vector fields Va, the flow of which uniquely specifles a 
congruence of curves along which the solution to the Gauss' law constraint 
can be constructed. The action of on any function Q{xi,Xj,Xk) is given 
by31 

^^Tho interpretation in terms of translation operators will notably facilitate the proof 
of convergence of the solutions, since for \Q\ < 00 Vf € E, one still has that l-ftTeQI < 
00 Vr G E for appropriate choice of Q. 

^^One must additionally have specified the functions corresponding to two eigenvalues 
Ai(r) and \2{r) of the CD J matrix '^ae, with the third eigenvalue Xsir) fixed by the 
Hamiltonian constraint. 

■^^ Assuming the convergence of the series (243), these two characterizations are equiv- 
alent. We assume convergence based on the Baker-Campbell-Hausdorff theorem, which 
guarantees the existence of an element C of the algebra, such that e'*e^ = e*^. 

^"Once a direction of integration has chosen and the boundary data specified, then the 
integration must be carried out without exception along that direction. 

^"This is due to the CDJ Ansatz, and holds for algebraically general spacetimes 

^^Not worrying about the covariant-contravariant position of the indices 
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^{Xi,Xj,Xk) = Va^Q{Xi,Xj,Xk) 



+exp 



<^{xi,Xj,0) 

0' 



- Pfkixi, Xj,0, Xj - d^ P%{xi, Xj, 0, Oj 

rx, qU - J^de /3-,{x^,x„e),x, - J.^de /3^(x„x„0 

f / — 7 7 7 

Jo (x, - /« d^' /3« (xi, xj,^'),Xj - /« d^' /3»,(xi, xj,^' 



5.2 Criteria for boundedness of the inversion 

It is not enough to prove convergence of the individual operators when 
inverting the Gauss' law constraint. One must additionally take into account 
the effect of the insertions Uae in (25), which also come into play in placing 
bounds upon the norm \^a\- This in turn depends upon choice of and 
the CD J matrix eigenvalues (Ai, A2), which are chosen Vx G S in advance of 
solving the constraint, and the source term Qa as well as the boundary data 
/3a- We prove well-behaved generic choices of Ai = Ai(r) and A2 = A2(r) 
and -B* = Bl^{r) place an upper bound on the solution. 

First, we require that for / = (Ai, A2, A3, S*) and their derivatives exist 
everywhere /(r) G C°°{T,) < 00 Vr G S. The following norms can be defined 
II / 11= max{|/(f)|} Vr G S and || /3 ||= max{|/3(r)|} We S. The zeroth 
order term of (25) is given by 

Go = /3(x, y,0) + K f = P + K T f{r). (248) 

Jo ^y{X, y ,Z) Jo 

The effect of the dual translation operator K does not affect the bound on 
I /I since it serves to translate the source within its domain of definition, 
which is chosen to be well-defined Vx G E based upon the configuration 
chosen. We require the boundary data to be similarly well-defined. Hence, 
one has that 



|Go| <|| ;S II + II / 



dy' =11 /? II + II / 



(249) 



The first order term is given by 
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^^There is no loss of generality by using the direction y, which is provided for clarity of 
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Gi = K dy — ^/3{x, 0,z) + K dy — — / —r- yj-^dy (.250) 

Jo Bl{x,y',z) Jo B^{x,y',z)Jo B^{x,y",z) 



Define the ratio 



"^''^ =c{f) (251) 



and the following norm || c ||= max{|c(r)|} Vr G S.^^ Then the second term 
of (250) is bounded by 



j\y'\c{y')\\Go{y')\ < J^' dy' \\ c \\\\ / || y' = |- || c |||| / || . (252) 
Then the following inequality holds 

<|| cllll /3 II 2/ + ^ II dill / II . (253) 
One finds the following recursion relation for the (n + 1)*'* term 

|G„+i| <|| P II + r dy'\c{y')\\Gniy% (254) 
Jo 

which leads to the following inequality 

II r II" 7/"- II f II /II ll"+l 

|G„+i| <|| ^ II + T^( ,V )• (255) 

n! II c II V (n + 1)! / 

Application of the Minkowski inequality to (25), where each succesive term 
is bounded by (255), one has 



II ^ ||n II f II °° II ^ ||"-+1 ,,"-+1 

(V + ar'Ql < V II /3 II + ^ E " ' .1 

^ c n (n + 1)! 

n=0 " " n=0 ^ ' 



=11 P II e"* + (e"""^ - 1) . (256) 



the exposition. \y\ = max{|a;|, \y\, \z\} simply plays the role of placing an upper bound on 
each term, which will in turn place an upper bound on the norm of (25). 

^^As with the previously defined norms, || c [| is also completely known in advance of 
carrying out the inversion based upon the specified configuration. 
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Hence, given the norms 



max 



max 



Q 



B 



V i,a Vr G S andS^ / 0; 
a,e,i,byre E andS* ^ 0; 
y[r\ = max{|x|, \y\, \z\} 



(257) 



where r = {x, y, z) in the last line of (257) corresponds to a particular point 
in S at which the constraint, is being solved, one has that 



|(v + a)"^Q| < oo (258) 

or that the full series (25) is bounded. Boundcdncss follows from the ob- 
servation that the norm of (25) is bounded by a pointwise convergent series 
(256). Note that boundedness for an infinite series does not imply conver- 
gence of the series. 
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6 Solution to the initial value constraints by ex- 
pansion about the Kodama state 



Let us now demonstrate a method to solve the constraints when there is a 
nonzero cosmological constant A. It is known that the Kodama state t/^Kod 
solves the constraints, which corresponds to a CD J matrix = ~x^ae- 
This automatically satisfies the diffeomorphism constraint 



eijkBiB^^ae = -jeijkBiB%e = 0, (259) 
due to antisymmetry of eijk- The Gauss' law constraint is given by 



We(*ae) = Bl^ + C/e(/a/,*,e + fefg^ag)- 



(260) 



For = -^6ae we have 



Bidi^ae + BiA{{Ufg^g, + fefg^ag) 
= BidiSae + BiA{{fafg6ge + fefgSag) 

= + BiA{{fafe + fefa)=0 (261) 

due to antisymmetry of the structure constants fabc- The Hamiltonian con- 
straint yields 



(Far* + Adet*) 



= 0. 



*ae=-(3/A)5a 



(262) 



The existence of an exact solution to the constraints motivates the search 
for other solutions by perturbative methods. So we will construct additional 
solutions by expansion about the solution "^ae = —j^^ae- Define the relation 



^ae = -{jSae + €ae), (263) 

where e^e parametrizes fluctuations about the conformally self-dual sector 
of Einstein manifolds^^ Substitution of (263) into the Gauss' law and dif- 
feomorphism constraints yields a cancellation of the —j^^ae term, and these 
constraints can be rewritten by replacement of with e^g to yield 



eijkBiB'^eaeO; We{e„e} = 0. (264) 
The indues and CP2 [?]. 
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The first equation of (264) implies that = e(ae) is symmetric in ae, whence 
Eae is now given by 



^ae = eLef + Ef,^f. (265) 

Note that the —^Sae term in (263) has no effect on the off-diagonal terms 
in (265). The Gauss' law constraint reduces to 



= Jj (266) 

where j/ is the Gauss' law propagator, which establishes a map from the 
anisotropy to the shear elements.^^ Upon substitution of (263) into the 
Hamiltonian constraint and cancelling off a factor of A one obtains the 
Hamiltonian constraint in terms of e = eae, given by 



A h? 

= tre + -Fare +—dete = 0. (267) 
o o 

Unlike for the kinematic constraints, there exist remnants of the cosmological 
constant A in (267) since the Hamiltonian constraint is nonlinear in ^ae- 
Observe that (267) involves only the invariants of eae, which can be written 
explicitly in terms of the eigenvalues (01, </>i, ^a) as 



2A A^ 

^l+<j)2 + h + — {4>l<f>2 + + 030l) + -^^/'102'?!'3 = 0. (268) 

One obtains a closed form solution for 03 in terms of (pi and 02, 



whence any reference to the Gauss' law constraint does not explicitly appear. 
Since the CD J matrix is symmetric according to (264) , when diagonalizable 
it can be written as a polar decomposition 



*6/ = -(e^'^)ba'(^<^a'e' +ei'e'01 +4e'02 +e^,,,03(01,02))(e ^■^)e'/(270) 



where 6 = [6^, 9"^, 6^) are a triple of complex 5*0(3, C) rotation parameters. 
While (270) encodes the Hamiltonian constraint solution within 0i and 02, 



^^Thus the shear elements are not independent degrees of freedom, but simply an image 
of the anisotropy elements. 
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regarded as the physical degrees of freedom, it does not incorporate the 
Gauss' law constraint Ga- In Paper VI it is shown how Ga can be used to 
fix the paramters 9 in the polar representation of the initial value constraints. 

On the other hand, when one attempts to solve the Hamiltonian con- 
straint in the Cartesian representation, there enter unphysical degrees of 
freedom. The Hamiltonian constraint (267) can be written as 

H = S'^eae + Al'^^'^^eaeebf + A2^»^''^'^^e„ee6/eeg (271) 
where we have defined the S0{3,C) invariant tensors 

jaebf _ l^^aejb/ _ ^eb^afy^ ^aebfcg _ ]_^abc^efg (272) 

on the 6-dimensional Euclidean space of of CDJ deviation matrix elements 
€ae- Upon implementation of the Gauss' law constraint one projects onto 
the 3-dimensional subspace of anisotropy elements eg using 

e„e = (eL + EUf)eg ^ Ti[A]ef. (273) 

One obtains a corresponding projection of the invariant tensors (272) onto 
corresponding vectors 

Qf'g' = jaebfr^frpg'^. Qf'g'h' ^ £;«e6/c3 _ ^274) 

Making the identification of the diagonal elements e// = e/, we can project 
the Hamiltonian constraint onto the reduced space as 

ei + e2 + e3 + Q[B„;ei,e2,e3] = 0, (275) 
which consists of a linearized part with a nonlinear error term Q, given by 

Qa = AQf^efEg + A^Qf^'^efegEh = QA[e; A]. (276) 

There are a few things to note regarding (275). First, it applies equally 
to the full theory as it does to the spatially homogeneous sector of A^. 
Additionally, it enables one to solve the constraints by fixed point iteration, 
by expansion in powers of the cosmological constant A for each configuration 
Af. Hence, one defines a recursion (e3)(n), where (e3)(o) = ~ei ~ ^2 is the 
zeroth order solution. Then one defines the recursion relation 
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(e3)(n+i) = -ei - £2 - QA[ei, £2, (e3)(n); (277) 
The full solution is given by €3 = lim„_^oo(e3)(n)) which can be written as 



*ae = - (jSae + T^e^i + T^.ea + Tle3[ei,e2; ^]) (278) 

Equations (278) signifies that the physical D.O.F. in the Cartesian repre- 
sentation are given by ei and 62, with the Gauss' law constraint encoded in 
the curved basis r/e.^^ Consistency of (278) with (270) implies the following 
relation 



Tie^f = {e'-^)aa'{e'-^)ee'ei,,,<Pf, (279) 

which relates the anisotropy elements directly to the eigenvalues. It is shown 
in Paper VI how (279) can be used to find the SO{3, C) angles 9. 



•^^Notc that the only part of the configuration variables in the homogeneous sector that 
enter into this expansion are ei, 62 and 63, which involve projections of {A~^)1 onto the 
direction of fe. The effect is to contract the spatial indices, which leaves remaining only in- 
ternal indices. This may be interpreted as fixing a gauge which breaks the diffeomorphism 
invariance, by singling out a special direction k in E. 
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7 Discussion 



In this paper we have placed the configuration space Ashtekar/instanton 
representation variables into a new context. It turns out in the latter vari- 
ables that the initial value constraints can be regarded as constraints on 
the momentum space of ^inst = i'^ae^^i)- When one rcsudcs this space 
to its physical degrees of freedom, this leaves remaining the question of the 
canonically conjugate variables on the full space. The association of the con- 
nection with a S0(3,C) frame and a choice of coordinates seems natural and 
it extends to the magnetic field i?*. On a more basic level we take the inte- 
gral curves 7 from whose vector fields these configuration variables can be 
derived as the fundamental quantities. Then each congruence determines 
an equivalence class of configurations on which the Gauss' law constraint 
may be evaluated. It then makes sense to regard the momentum space 
variables on the constraint shell as being aspects of the instanton represen- 
tation which are independent of coordinates and of the S0(3,C) frame. We 
have performed a prc-ciirsory analysed various configurations based upon 
the integr ability of the vector fields, developing intuition on certain avail- 
able configurations. Once one finds an integrable configuration, there seems 
to be no obstruction to pairing it with different values for the CDJ matrix. 
We have found a large class of configurations for which the Gauss' law con- 
straint is integrable. This provides a large class of configurations for which 
one can in principle construct solutions to the Einstein equations according 
to the prescription of Paper II. 
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8 Appendix A: Expansion of the helicity density 
matrix in the tensor representation 

We would like to write the Gauss' law constraint in a form more amenable 
to solution, in conjunction with the Hamiltonian and the difFeomorphism 
constraints, to the initial value problem of general relativity. Rewriting the 
constraint for completeness, we have 

We(^'ae) ^ ^I'i^fg) = S^^M^fg) + C^^^fg = GQa- (280) 

To compute the matrix elements of the connection C/^ we will need to 
evaluate the components of the matrix product C^e = BIA'- in the tensor 
representation.^'^ The SU{2) (g) S'?7(2)-valued matrix Cbe constitutes the 
connection corresponding to the S'J7(2)_ Gauss' law covariant derivative 
of the CDJ deviation matrix elements. The tensor representation of the 
connection for this covariant derivative is given by 

{Cays = Cbe {fabfSge + febgSaf) (281) 

where the structure constants are given by fabc = ^abc- We will now explicitly 
compute the components. 

Starting with the 'even' elements corresponding to the first component 
of the SU{2)_ charge 

(Ci)^^ = Cbe{eiblS2e + eeb2Sll) = C13 — C31 = — C[3i]; 
(Cl)^^ = Cf,e(eife253e + eeb3<5l2) = ^336132 = — C33; 

(Ci)^^ = C'fte(ei635ie + ^ebl^is) = ei23C'21 = C'21- (282) 

Moving on the the 'odd' elements corresponding to the first component of 
the SU{2)^ charge, 

(Ci)^^ = (ei62<5le + ee61'^12) = £132(^31 = — C31; 
(C^)32 = Cfee (ei63(52e + eeb2Sl3) = €123^22 = C22; 
(Ci)^^ = Cbe{€lblS3e + ee63'^ll) = C'21 — C12 = — C[i2]. (283) 

Moving on to the 'diagonal' elements corresponding to the first component 
of the SU{2)^ charge, 

Observe the order of the factors in this product, as well as the convention for lowering 
the indices. 



72 



{ClY^ — Cbe{eiblSle + Cefcl^ii) = C32 — C23 = — Cps]; 
(C'l)^^ = C'6e(ei62^2e + ^eb2Sl2) = C'32ei32 = —(^32; 

{Cif^ = Cbe (eiba^ae + eebsSis) = Csaeisa = C23. (284) 

Moving on to the 'even' elements corresponding to the second component 
of the SU{2)- charge, 

((^2)"^^ = C'6e(e261'^2e + ee62<^2l) = £231(^32 = (^32; 
((^2)^^ = C'be (e262'^3e + ee63<^22) = ^"21 — (^12 = —C[i2]', 

{C2f^ = Cbe{e2l,3Sie + Ceblfe) = e213Cll = -Cn- (285) 

Moving on to the 'odd' elements corresponding to the second component of 
the SU{2)^ charge, 

((^2)^^ = C'6e(e262'^le + Ceblfe) = (^32 — (^23 = —C[2Z\\ 
((^2)^^ = Cbe(e263<^2e + ee&2'^23) = 6213(^12 = — C12; 

((^2)^^ = (^66(^261 <^3e + ee&3'^21 = €231(^33 = (^33- (286) 

Moving on to the 'diagonal' elements corresponding to the second component 
of the SU(2)^ charge, 

(6*2)"^^ = Cfee(e261'^le + eefel'^2l) = £231^^31 = C'si! 
((^2)^^ = Cbe{f2b2^2e + £662^22)^13 — C31 = — Cpij; 
((^2)^^ = C'fte (e263(^3e + ee63<^23) = C'l3e213 = — (^13- (287) 

And finally, we move on to the elements corresponding to the third compo- 
nent of the SU{2)- charge. Starting with the 'even' components, we have 

(^3)"^^ = Cbe{^3bl^2e + ee62'^3l) = C'22e321 = — ^^22; 
(Cs)^^ = C'6e(e3fe2<^3e + ee63<^32) = ^SuCls = C13; 
(Cs)^^ = Cbe{e3b3^le + ^eblhs) = C32 - C23 = -C[23]- (288) 

Moving on to the 'antisymmetric' elements corresponding to the third com- 
ponent of the 5J7(2)_ charge, 

(C's)^^ = Cbe{e3b2^1e + ^ebl^32 = C'lie312 = Cu] 
(Ca)^^ = Cbe{^3b3he + ^eb2^33) = C13 — C31 = — C'^ij; 

(Ca)^^' = Cbe{e3blS3e + eeb3S3l) = C23e321 = -C23. (289) 
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And finally, moving on to the 'diagonal' elements corresponding to the third 
component of the SU{2)- charge, 



{C2 



\33 



(C3 

= Cbe{^3b3^3e + ee63'^33) 



C21 — C 



12 



- C12; 
-C[12]. 



(290) 
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9 Appendix B: Results from the standard config- 
urations 



9.1 Abelian contribution 

We now compute the constituents for the Gauss' law constraint for a wide 
class of physically interesting configurations. We will display the intermedi- 
ate steps of the calculation to eluciate on how various structures of interest 
recur. Take the Ashtekar magnetic field to be of the form^^ 



(291) 



where a = a{t), /3 = j5{t) and 7 = ^{t) are spatially constant. First we will 
compute the Ashtekar magnetic field for (291), given by 



c 



(292) 



Starting with the terms involving spatial gradients we have 



^ijkdjAl — dj 



a 



(3 



(OiX"^ \ / / 1 \ / X^X"^ \ 

+^[-^y^'^' + y,{^a,x'^ + ^k,x'^) 



(293) 



This simplifies to 



(294) 



In the first line of (294) the second them cancels with the second member 
of the third term in brackets and the fourth and sixth terms, in the second 
line, vanish due to antisymmetry of e^^^ . The result is that the gradient 
contribution to B^, the 'Abelian' part, is given by 



(295) 



•''**Wc take the connection to have mass dimension of [Af] = 1, as expected for a gauge 
field. This fixes the position dependence of the individual terms. 
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9.2 Nonabelian contribution 

The nonabelian contribution to is given by 



2 ^ ^~ 2 



a 



,3 



X 



Equation (296) splits into nine terms. The first term is given by 



(296) 



^e'^''r^%kSbj (^) = Sai (^) . 



Moving on to the second term of (296) we have 



(297) 



-eiamx-(|§)(.298) 



We have used the vanishing of the first term on the right hand side. Moving 
on the third term of (296), we have 



1 ijfc 
2 



V2r4/ 



(299) 



where we have used = r^. 

The fourth term of (296) is given by 



le^'T'^e^knSbjX-l^) = ie,,,(<5,"<5;;-<5^5^)<5,,x"(^) =6-v(^)(.300) 
Moving on to the fifth term of (296), we have 



= ^(«-«-M + 3«)x-x-(^) =a;«x'(^). (301) 



1 

_ijkfabc T.nrf.mr_ 
J ^ckn^bjm-^ •'^ Xf'^ 



Moving on to the sixth term of (296), we have 
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± ijk fabc b j(JP\ 

2 ^ J ^ckn-'^ X J. V / 

= - K^Dx^^'x^ (^) = ^^'^x^ (^) . (302) 

The seventh term of (296) is given by 



Moving on to the eight term of (296), we have 



\jjk fabc . m c k ( 1^1 \ 

= l^^^'m - ^?'m)x-x'^x''(^) = -e''"'x''{^) (304) 

where we have use the fact that the first term vanishes. Finally, moving on 
to the last term of (296) we have 



l^ijfcy-abc^b^j^c^fcj^^!^ ^0. (305) 

Combining the results of (295) and (297) through (305), we have that the 
Ashtekar magnetic field for (291) is given by 



B: = 6-^^^^ + xV[^^^^) - e-V^^^. (306) 
The vector fields = B^di for the configuration (291) are given by 



where d/dt = Xidi and Ja are the angular momentum operators. 



9.3 Magnetic helicity density matrix 

The magnetic helicity density matrix Cae, which can be thought of as the 
direction cosines between the the coordinate axes of A'^ and , is given by 
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Cae = AfBi = [6ai (^) + CaijX^ (^) + x"^ 



X 



./3(a + 7) 



Equation (308) splits into nine terms, given by 



(308) 



^^.^e «7(« + 7) ^ ,a,e 7(-a7 + /3^) _ o_ 



(309) 



Collecting like terms, this reduces to 



Cae = S" 



+■ 



^ ((« + /3)(/32 + a^) + (-«7 + /32)(« + 

' af5{—a — /3 + a + 7) 



) 



^ ^ae (Q^-/3^)(a+7) ^ g^a^e ^^(« + 7) ^ ^aem^m «/3(7 - Q) 



(310) 
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10 Appendix C: Roots of the cubic polynomial in 
closed form 



We would like to solve the cubic equation 



+ Ae'^ + Be + C = (311) 

for some coefHcients A, B and C. First we perform the transformation 
e = y + r, yielding 

(y + rf + A{y + rf + B{y + r) + C 
= y^ + (3r + A)y'^ + (Br^ + 2Ar + B)y + + Ar'^ + Br + C = 0. (312) 

By the method of Cardano we choose r = — ^ in order to eliminate the 
quadratic term in (312), which yields 



y +py = q 



where 



(313) 



V 3 27 / 
making the further substitution y = (4p/3)^/^m, we have 



(314) 



m + -m 
4 



4p\-3/2 



(315) 



Letting m = sinO and making use of the identity 4sm^9 + SsinO = sinSO, we 
have 



sin3^ = g(^) ; 6 = \sm-^{4q{4p/3)-^^^ + p) 



(316) 



.4px -3/2 

where p = ^ labels the the three roots for n = (0, ±1). So the solution of 
(315) is given by 



m = sinO = T^^g 



4g(4p/3)-^/2 



(317) 



where we have defined Ta{x) = sin[asin ■^x + p]. This leads to the general 
solution to (311) of 



A /4/„ A3\^^ 



-4(C 



AB I 2A^ ' 
3 27 - 



3/2 



(318) 
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